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ABSTRACT 
A new method for design of control systems in the time domain is 
presented. The adjustable parameters, the roots and the zeros of the 
system are evaluated simultaneously either by hand or by computer. 
The method is applicable to any system configuration and is able to 
handle a large number of adjustable parameters. The method is analy- 


tical, any degree of accuracy can be obtained. 
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PART {I ANALYSIS 
CHAPTER I 
THEORETICAL CONSIDERATION OF ALGEBRAIC EQUATIONS 
1-1. Introduction. 

The response of feedback control systems is governed by the complete 
set of the closed loop pole-zero configurations. Since the zeros are 
usually known, the only problem is the poles. For linear control systems, 
the problem of poles is conveyed to the roots of the characteristic equa- 
tion which is an algebraic equation. Several methods have been developed 
to this end such as Guilemin's mached:. root-locus method and Mitrovic's 


Mernod- >. 


Guillemin's method starts the synthesis by choosing a set 
of pole-zero configurations from the consideration of specifications and 
then finds out the network needed to give this chosen pole-zero configura- 
tion. The root-locus method is essentially a graphic method for one vari- 
able which has linear relation with the coefficients of the characteristic 
equation. Mitrovic's method is also a graphic method but handles two 
variable coefficients, The method proposed in this paper is rather an 
analytic method and is able to handle any number of variables which are 
contained in the coefficients of the characteristic equation. The analy- 
sis is based on the root-coefficient relation of an algebraic equation, but 
a new technique of partitioning the equation is employed. 
1-2. Root-coefficient Relationship. 

For 4 ae order algebraic equation, there are n relations among the 
roots and the coefficients of the equation. If all the coefficients are 


actually known numerically, then all the roots are fixed in such a way that 





they satisfy the root-coefficient relations. Conversely, if a set of 
roots of the characteristic equation are chosen, then all the coefficients 
of the equation are fixed by those relations. However, if variables are 
contained in the coefficients, then the roots of the equation are varied 
in a definite relation to the variation of the variables. For example, 


if one coefficient of the characteristic equation is variable, then the 





relation between the roots and the variable coefficient can be plotted by 
root-locus method. But if this relationship is locked at from a different 
point of view, an analytical relation can be obtained. Consider a third 
order characteristic equation with one variable coefficient as shown in 


the following equation: 


S°+ BS + f£(4)S +B, =O (1-1) 


where Bs and B, are known numerically, f(A) is a function of % 
which can be varied. Assume — -S, and —53 are the three roots of 


the equation, then the root-coefficient relations are as follows: 





Se ae en (1-2-1) 
SS oS Sey mca Taz) 
esr or ai). (1-2-3) 


Since f(A) is adjustable, then equation (1-2-2) can always be satisfied 
for any choice of the rocts by varying filb) accordingly. However, 
equations (1-2-1) and (1-2-3) set constraints to the value of the roots. 


Namely, the three roots must satisfy the two equations (1-2-1) and (1-2-3). 





If the three roots are treated as the variables of a set of simultaneous 
equations, then there are three unknowns but only two equations, consequent- 
ly one root is arbitrary. This arbitrary root can be chosen independently 
if there is no restriction to the value of t (fh) , and therefore this 
root can be treated as an independent variable in the root-coefficient 
relations. In the root-locus method, the independent variable is f.(f) 
while in this point of view, the independent variable is conveyed to one 
of the roots of the characteristic equation. 

The above reasoning is based cn the solution of simultaneous alge- 
braic equations, i.e., Sy S, and = are treated as independent 
variables. Since those variables are the roots of an algebraic equation, 
a complex root must be treated as two variables because of its conjugate 
root. If one variable of a pair of complex conjugated roots is expected 
to be treated as the independent variable, this variable cannot be chosen 
to be the root itself, but can be treated as follows. Assume oe and = 55 
are a pair of complex conjugate roots of the equation (1-1), those two 
roots can be transformed to other two variables by the following trans- 


formations: 


2) SS ae 2PWny (1-3) 


C 
5,52 — On (1-4) 


Where 5 and Wry are known as the damping ratio and the natural frequency 


of the pair of complex conjugate roots —5, and -S, . By this trans- 


2. 

formation, both g and UJ), can be treated as independent variables. 
In equation (1-1), if f(A) was also a fixed number, then all 

roots were fixed numbers by equations (1-2) and no root was arbitrary. 


The arbitrary roots arise from the introduction of variables to the coeffi- 


cients, If one variable is introduced to the coefficients, it has been 
3 





shown one root is arbitrary. For the same reason, if two variables are 
introduced to the coefficient of a ia order characteristic equation, 
then two roots are arbitrary, and for MM variables to the coefficients, 
IM xroots are arbitrary. The root-locus method is the case of one vari- 
able while Mitrovic's method is the case of two variables, but both methods 
are quite different ways of approach. [It is because of this way of look- 
ing at the root-coefficient relations, an analytic method can be derived 
and is capable of handling any number of variables contained in the co- 
efficients. 
1-3, Partition of a Characteristic Equation. 

If adjustable parameters are introduced to the coefficients of a 
characteristic equation, then this equation can be partitioned into two 
equations. One of them is formed by the arbitrary roots, the other is 
formed by the constrained roots. The constrained roots here are defined 
as those roots other than arbitrary. 


Consider a 5th order characteristic equation 
& a 3 2 
Seen, tees + Bee 5, oe OC Sen!) 


Assume —5, -S,  -5, - 5, and -S, are the roots. The negative 
, J J 


Sign is for convenience. By expanding the equation 
eect cst >. S452) (St5e) == © (1-5-2) 


Set the corresponding coefficients of (1-5-1) and (1-5-2) equal, then the 


root-coefficient relations are obtained as follows: 


B, == S$, +52+S3 +5445 Green) 


pe = ie.) (53 Tet, Fes. +(5,t-, tO, oe (1-6-2) 








~~ 


Bee 8.5,(5,tS2) + (5, 44,) 5,5. t ere eS.) 4 Ss, SS. (1-6-3) 


Zz 


o,f 1 (S,45.) 55, + 5,5. (5, +5,) | Se (1-6-4) 


Bo = 3,52 5354 5s (1-6-5) 


Assume the coefficients B, and an are adjustable, then equation (1-5-1) 


can be written as follows: 


ees + BS 3 (4, AL)S + FG. f)s4+B,=0 (-7) 


Where 7 Gi. hi.) and Teal f, . h, ) are functions of a and 2 which 
are adjustable parameters. By : B and Be are assumed given fixed 
values. A terminology, "controlled characteristic equation" is given to 
equation (1-7) on the reason that its roots are controllable by the ad- 
justable parameters. Since there are two adjustable parameters intro- 
duced to the coefficients, two roots can be chosen arbitrarily. Let those 
two arbitrarily chosen roots be ( -S, ) and (-S,), then the other three 
mors (->,), ¢ oy and ( ~DS— ) are not arbitrary but constrained to the 
fixed coefficients ( ee ; B, and a ) and the arbitrarily chosen roots 
as shown in the last section. By separating the arbitrary roots and the 
constrained roots, equations (1-6-1), (1-6-2) and (1-6-5) can be arranged 
as follows: 
From (1-6-1) 
S,+ 4,45, = (GS as Sy, (1-8) 

From equation (1-6-2) 

B, = ( 5, + S, +55)(5,452) + (5,544 S5¢# Se 52) +55, (1-9) 
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From equation (1-6-5) 

Bo = (5,5.)(7;%.55) (1-10) 
It can be recognized that ( = + 5, + S¢ ea Ee oe t Se S; ) and 
( S,S,S5, ) are the coefficients of a third order algebraic equation 


which has roots of -5, » —y, and -5S_. Define 


> 
Ce 5 >; (1-11-1) 
ees. + SESE (1-11-2) 
Boi Se Sy t Se Ean) 


Then the equation of the constrained roots is formed by 


s+ C,s°+ ise — © (1) 


It must be noticed that equations (1-12) has the roots -5,, ~S, and aoe 
Substitute the definitions of (1-11) into (1-8), (1-9) and (1-10) and 


manipulating, one obtains 


Co se C5 7 Ss) (eee) 


1-12-2 
e a Pete Smut) S,S. 


Gee = Gb: Vie (1-12-3) 
Consider equations (1-12), C,, (C, and Cy, are functions of the arbitrary 
roots and the fixed coefficients 6B , 68, and Ba of the controlled 
characteristic equation (1-7). Therefore the roots of equation (1-12) 
are also functions of the arbitrary roots ~-S, and =o and the fixed co- 


efficients. Assume there are no restrictions to the variable coefficients 
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it, f,,) and f. (hh , fh, ) then the controlled characteristic equation (1-7) 

is partitioned into the following equation 
(S+5,) (Ss eee tC, el Sioa &) li 

because both equations have the same roots. As Cs a, and Cc. are 
functions of S and S, , then the five roots of the original character- 
istic equation (1-7) are completely determined by equation (1-12). Equa- 
tion (1-12) is a third order algebraic equation and the coefficients are 
expressed analytically by (1-12-1), (1-12-2) and (1-12-3). Then the 
original 5th order characteristic equation is now reduced to a third order 
equation and by this reason, a terminology ''REDUCED CHARACTERISTIC EQUATIONS 
is defined for equation (1-12). 

In the above arguments it is assumed that the variables h, and he 
are adjusted such that two arbitrary roots are chosen. Therefore, the 
variables h and h. , which are the actual physical adjustable parameters 
of the system, must be determined. From equation (1-6-3), and (1-6-4) 
replace 6, by ie (h,, h, ) and 6, by cf ( hy, h,) since they are 


assumed adjustable, and manipulating, one obtains 


sak h, ; fz) ia C. + G, ( — +S2) fe G S,52 (Tee) 


f, Gh, fa) ==. 0.05.4 42) +O, S52 Cele) 


Substitute ee » C, and GS from equations (1-12) into (1-14), obtain 
= 3 1-15-1 
f.(h,,8,) = 8/55, + (545) - Be (5,15,) + 552 tS 45:)-255, (545)! 


2 
ee ae) == B05 tS) /s,S. + B,S,5. — B,5,5.(5 t5,) + em 5/5.) leS2) 


It can be seen from equations (1-15) that i, (4, k,) and fF, (h,, 43 ) are 


expressed analytically as a function of the arbitrary roots and the fixed 








coefficients ee , Bs and By, - And 4, f,also can be expressed ex- 
plicitly as functions of the arbitrary roots. Therefore if a set of 
roots has been chosen from the consideration of the reduced characteristic 
equation, then the physical adjustable parameters are readily determined 
by equations (1-15) or (1-14). 

By this partitioning process it can be seen easily that the original 
fice root-coefficient relations of (1-6) are transformed into other five 
relations which are equations (1-12) and (1-14) or (1-13) and (1-15). 


For convenience (1-12) and (1-14) are repeated here: 

B, eo Crs ct 5) _ 5,S2 

‘ B. sks (1-16) 
Fo = Cot CS tS) +O SS 


| 


| 


6 
C. 
& 


\ 


te (Cie (oS, te) + Gi Ses, 


In actual numerical computations, it is simpler to use equation (1-16) 
instead of (1-12) and (1-15), because of the progressive nature of the 
root-coefficient relations of algebraic equations, Namely, compute (C, 
and C, first, then using these calculated data to compute C, and then 
af and fos : 

In general, the n root-coefficient relations of a nth order character- 
istic equation which has variables in coefficients can be arranged in any 
form as one wishes. The process derived here has certain significances 
in the design of control systems. First of all, the order of the character- 
istic equations which determines all the roots of the system is reduced by 


the number of the variable parameters. Secondary, the physical adjustable 





parameters are expressed directly by the arbitrarily chosen roots. The 
stability and response of the system is determined completely by the re- 
duced characteristic equation while the adjustable parameters are deter- 
mined directly by the arbitrary roots. Since equations (1-16) are alge- 
braic equations and if no singularity is involved, then solutions of by, 
and h which actually determine the physical adjustable parameters al- 
ways exist. Although this process does not guarantee the adjustable para- 
meter obtained in this procedure to be always realizable, the realizable 
region can be found. In the design of a control system, this is a logical 
procedure, namely the pole-zero configuration of the closed loop system is 
chosen first from the specifications, and then compute compensator para- 
meters directly from the chosen pole-zero configurations. For every 
choice of pole-zero configuration, a corresponding set of values of the 
compensator are obtained. This partitioning process effectively trans- 
forms the desired pole-zero configuration to be the independent variables 
of the design procedure. Moreover, most of the designs are based on a 
pair of dominant complex conjugate roots, then two of the arbitrary roots 
can be chosen (not necessary) as the dominant roots. The adjustable para- 
meters of the compensator therefore can be expressed explicitly as func- 
tions of the dominant roots if the dominant roots are defined. 

In this process, there are no limitations. For a control system of 
any order, any number of adjustable parameters and any type of compensa- 
tor; the same technique can be applied. 

1-4. General Procedures of Partitioning a Controlled Cheracteristic 

Equation. 

The example in the last section assumed that the adjustable parameters 


are contained in two coefficients. Now consider a cascade single section 





R-C compensator as shown in the block diagram of Fig. 1-1, in which the 
compensator parameters enter all the coefficients. The characteristic 


equation of the compensated system is 


st4 (Btp)S'+ (8+ B ps +(BptxBys+PB=o0 (1-17) 
Mee b=K, B=-ab, B.=—atb 


j 


O alee 2) K 
bs (Sia S(Sta)(stb) 


Fig. 1-1 Single section cascade compensator. 
x= P/z ,; Ks a and b are know numerically. 
Assume -5 , Se ~S3 and 7 2 are the roots of (1-17). By expanding 
the equation 


oe ta) Ss) =) = © (1-18) 


and setting the corresponding coefficients of (1-17) and (1-18) equal, 


one obtains the following four relations: 


Bot Pao Set S, hee (1-19-1) 
EP imoeiae Get 2) Goat on at So. 6.5.5, (1-19-2) 
B,P+ 4B, = 5,52( 5; +54) + $3%(5,45 ) (1-19-3) 
PB, = 3,71 53 Sz. (1-19-4) 


Since there are assumed two adjustable parameters X and fe of the compensa- 


tor, there are two arbitrary roots. Assume -S , and oS are the 
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arbitrary roots, then ( en 24 ) and ( Doe. ) are the coefficients 


of the reduced characteristic equation. Define 


——. (1-20-1) 


ms | 
Cee ee (1-20-2) 
then the reduced characteristic equation is 
2. aes 
Se cy> = Co OC (1-21) 


Substitute the definitions of (1-20) into (1-19), one obtains 


| 


— | _ 

| O ne peat) (1-22) 
ee (515) ——2, —- 

| O ae a B O 


O 
| (5,45) | -B, OG — 5,2, 
B 
| 
. ie 


»| | p O 
x 

Equations (1-22) are the partitioned four root-coefficients relations in 
which C, and C. define the reduced characteristic equation while [Pp and < 
are the compensator parameters. These four dependent variables (Cys Cos 

Dp and & ) can be solved in terms of the arbitrary roots and the fixed 
plant parameters ( B, » B, and B,). For numerical computation purposes, 
only one dependent variable need be expressed explicitly as a function of 


the arbitrary roots, the others are computed readily by the previous cal- 


culations. Here C. is solved from (1-22) as shown in the following equa- 





tion 
| 2 
Corn | 3,2, —(5,+5,) —B, + B, (5, +5.) | (1-23-1) 
where 
Le | B,—(5,+52) | = — {| (1-23-2) 
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Now consider another example which introduces three adjustable 
parameters, Fig. 1-2 is a fourth order system with first and second 


derivative feedback. The parameters hy» hy and K (= KK) are assumed 





Ke 





Fig. 1-2 Pure derivative feedback K, h 


1 


and h, are variables where 


K = KK 


adjustable. The characteristic equation is 
s'+ BS + (B+ Kh)s4+(B,¢kKh)St+tkK =O (1-24) 
where B,= atb+c , Bee Ob Phere , and B =abdc ; 


Assume the roots are -S, , -S,, -S3 and - 5, then the root coefficients 


relations are as follows: 


a 5 Poros S, (251) 
Eien (Sah Se Sy SiSaF (5, #505, ) Se (1-25-2) 
eye Kye 5,5, StS, | 5,5, + $.(5,4+5,)] (1-25-3) 
KK = wy) S, Sontg (1-25-4) 
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Since there are three adjustable parameters (h, hy and K), three roots 
are arbitrary. Let these arbitrary roots be =) : - 5, » and - 5, : 
The reduced characteristic equation is of order one and its roots ieee 


Define Sy = Cos then the reduced characteristic equation is 


See =e (1-26) 


Substitute the definition of Sy = Co into (1-25), the root-coefficients 


relations become 


Gee ee eo, to tS, ) (127-7) 
Kon, == (3,4 5,)S5;, + S,S, + Co (5,45, +5,) — Bz (1-27-2) 
eee Ss | OS + ( 5, + 55) S01 @lecgeey 
K = >%29;Co (1-27-4) 


Equations (1-27) are the partitioned root-coefficient relations in which 
(1-27-1) determine the system response while the other three equations 
determine the adjustable parameters so that to give the roots determined 
by (1-27-1). 

From the three examples illustrated above, it can be seen that this 
partitioning process can be applied to systems of any order, any number of 
adjustable parameters and any system configurations. The important concept 
of this process is to separate the roots of the controlled characteristic 
equation into two parts: The arbitrary roots and the constrained roots. 
In the n root-coefficient relations, the arbitrary roots are treated as 
independent variables while the constrained roots and the adjustable para- 


meters are the dependent variables. The constrained roots are obtained 
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from the solution of the reduced characteristic equation while the ad- 
justable parameters are calculated from the arbitrary roots directly. 
The general procedures are summarized as follows; 

(1) Formulate the controlled characteristic equation’. Set 
the coefficients equal to the corresponding root-ccefficient. 

(2) Define the arbitrary roots+ and the coefficients of the 


reduced characteristic equation (Cos Creee)e Substitute Cos C into 


ve’ 
the equations in (1). 
(3) Solve for the aijustable parameters and C's in terms of the 

arbitrary roots. 
1-5. Summary. ; 

This partitioning process effectively transforms the independent vari- 
ables from the adjustable parameters to the arbitrary roots of the character 
istic equation. As the adjustable parameters (compensator) are determined 
from the desired closed loop pole-zero configuration, a design procedure 
for control systems based on the time domain can be derived. The parti- | 
tioning procedures presented in the last section are applicable to all 


system configurations. For specific cases of a nth order system, general for- 


mulas can be derived. For the cases that some of the coefficients themselves 


*In discontinuous operation (2), the controlled characteristic equations 
may be different for different modes. As far as this step is concerned, 
each mode must be treated independently. 


+ The number of arbitrary roots is equal to the number of the independent 

adjustable parameters. The irdependent adjustable parameters may be less 

than the compensator parameters because of the specifications and other 
restrictions. For details, it is discussed in Chapter IV. 
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are treated as variables, general expressions from the partitioning 
process have been derived as shown in Table 1-1 to Table 1-6. For 
other cases, they are shown in the subsequent Chapters. In order to 
illustrate the applicability, the following examples are presented. 


Example: Given a characteristic equation 


5s? + 1oS*+ j285°+ 520S8° + 13008 + 3000 =O (1-28) 


which is unstable. Assume a network can be devised such that the co- 
efficients of the first, second and third derivative terms can be ad- 


justed. Then the controlled characteristic equation is 
4. 3 ~ | os 
Peerios + 7.5 + 7,5 + f,5 + 3000=0 (1-29) 


where an ; he and 7 are variables. From Table 1-5, the coefficients 


of the reduced characteristic equation are: 


Geen ( 545,75.) (esol) 


Co = Bo /5, ee (1-30-2) 


The functions of the variable coefficients are: 


ome C( So, + 5,5, + 5:5,) +C, 5,5 S2 (230-5) 
ie eG 5, FS) S,) + C, (5,5. 45.5, + 5:5,) + 5.5253 (1-30-4) 
1 = Co iy eSjceeds fS5) ef (Sak o> i Sy Sy) (1-30-5) 


where (-S, ), ( -5, ) and ( - 5; ) are the arbitrarily chosen roots. 


Since equation (1-29) is of the 5th order, one root must be real. Assume 
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this real root is C2 For convenience let (-5S, ) and (-S,) be 


the dominant complex roots and express them by the following notations: 


S,+5, = 25 Wn 


zi 
So 


(1-31-1) 


ie ad=2) 


Substitute (1-31-1), (1-31-2), By = 16, and Bo= 3000 into the equa- 


tions (1-30), one obtains: 


ee or (2, ts) 


aa 3000 /Wn S3 


| 


f= Co( Wertz $wo5s) + C\Wn Ss 
1s = Ge 27W, t Sy ale C, (w,, + 2 Pun 5s) cf Ww, Ss 


fe == Co tG,(2fwWat Ss) t (Wr +29 Uns ) 


The reduced characteristic equation is 


i C,5 1G5— © 


(1-32-1) 


(123222) 


(1-32-3) 


@le32=4) 


(=3255)) 


less) 


which determines all the roots of equation (1-29) and consequently deter- 


mines the stability boundary. By putting all the three arbitrary roots 


on the left-half plane, then from (1-37 2) Co is always positive, the 


condition of stability for the arbitrary roots therefore from (1-32) is 


C,>0 


ae llo = 2PWy + 3 


16 


(1-33) 





The inequality (1-33) implies that the three arbitrary chosen roots are 
confined in a certain region if the system is to be stable. This region 
is plotted in Fig. 1-3. Notice GWy is the real part of the complex 
roots, If S, and S) are complex the boundary is a line ( S = -8). If 


1 
Sy and S5 are real, the boundary is a point, (S = -8 + jO ye 
Assume the dominant roots (-5S, ) and ( ~5, ) are desired such 
that all other roots must be to the left of the dominant roots, then 
the following inequalities must be satisfied 


>, = FW (1-34-1) 


and een (1-34-2) 


Substitute the boundary conditions of (1-34-1) and (1-34-2) into (1-32-1) 


one obtains 


On = 16/5 (1-35) 


Equation (1-35) defines the dominant root region which is also shown in 


Fig. 1-3. Within this region, choose i = 0.5, Wh= 4, S43 = 4, then 
from (1-32-1) and (1-32-2) 

Cc,= 8 

Co — 47 
The roots of (1-33) are: “44 ) 5.6 ( W, = 6.85, § = 0.59) 


The root configuration for this choice of the arbitrary roots is also 
shown in Fig. 1-3. 
The next step of the design is to calculate the variable coefficients. 


Substitute § = 0.5, W), =4, S, = 4, C, = 8 and C, = 47 into equations 


3 i 
(1-32-3), (1-32-4) and (1-32-5), one obtains 


+, = 2015 
tT, = 696 
f3 = 143 
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For another choice of the arbitrary roots, the coefficients of the 
reduced characteristic equation and the variable coefficients are com- 
puted in the same way. 

The above example assumed that three coefficients are variable. Other 
cases, such as a cascade compensator, feedback compensator, etc., can be 
carried out in the same procedure, since the expressions of the coeffici- 
ents of the reduced characteristic equation and the adjustable parameter 
are obtained readily. The stable root region an*dominant oot region are 
important in design, the general methods for evaluating those two regions 


are discussed in the next chapter. 


1g 
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F, = Co (5,452) tO, 5S. 


f,= Co+C,(5,45,)4GS,S, 






— ee ($,+5.) 




















th C, = B-S (5,45) -S.% 
C= B, ~ C, (5, +5,) a $152 
— Be 
Co ie $,S2 , 
Can= | 
Cn-3 — Bp, —C,, (s +5.) 
Cy-4 = Baz —Ch; (5,45) - G25, 
NM th Co-t = Bn3 —Crrg (545)-Cn355, 






en= B,- Cz (5,t5,) -C3z 5,52 






Table 1-1 Functions of the variable coefficients f, and fz . The. 
functions contained in this table are assumed only B, and B, are adjustable. 


(-S,) and (-S,) are.the arbitrarily, chosen roots and'C's aré the coefficients 
of the Podueed characteristic equation. | 


s+ ais” 4 at Jee. — O 


7a; 
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C, = [6,— Co(5,+5,)] AS 
Ce — i. /'3,5; 

















e — 6, ( 5 + s,) 
C, = [8 —-065,45,)] 432 
Ca Bo Ve SER 


Red 

s — Bes Cs fsa) 

e. = By - Cz (S5,¢ 5.) - S52 

ee | BE ICICS +5.) sion 






$,= PCs) +G 55s 






jee G,(5,t5,) 10,55. 









Cy-2 = | 
C3 — Bae a ¢ t5:) 
Cy4= Baz — Cn-3 (5,t5)-Cry 55a 






ie 
| 


B,— C; (5,+5,) —Cy 5,52 
[ B.~ C0t5,45,)]/ S52 
Bo / 3:52 
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Table 1-2, Functions of the variable coefficients i wand fs; 
The functions contained in this table are assumed only B, and B. are 
adjustable. (-S_) and (-S,) are the arbitrarily chosen foots and C's 
are the coefficients of the reduced characteristic equation, 
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_ C=, 
ce — ( Be Gc -5,S,) Ast52) 


“F th 
= SE Sm, 







C, = B,—-(5,+ 52) 




















Sth 
= ( eS Oss ce 
— B, /'S.S, 
t= Co(5S,+5,) + GS 52 
Be (Sioz) t= C, + G(545,) +C; $52 
© th GC, = By —Cs (5,45,)— 5,52 _ 


oe, (8 -C, — G 518.) /S45;) 
a. / 5,52 












Ces = B. = Cp-2 ( 5,45.) 
~. Cn-4 — Ba aa Cn; ( 5,45.) ~€n-25,52 















B, -C, (5,48,)-GQ 55 
[ B-G@—-G55.] 55.) 
Be /S,S2 






HW tl 






Table 1-3 Functions of the variable coefficients fi and f;. 
The functions contained in this table are assumed only B, and B, are 
adjustable. (-S,) and (-S,) are the arbitrarily chosen foots and C's 
are coefficients of the reduced characteristic equation. 





order coeff. of the reduced Functions of the 
ae chara. equation variable coeff. 
equa 

3rd. 


Sore Gi) t 22) Gans, Ss 
4- th C= B, — C2 (5,+ 52) 





C= | 
Ga Ce (5 +5,)C aS See te Co 5, See 
C, —= B; aes tS. ) Cy ae 5,52C3 
oe pres) fapserco. +5, ) FEn5 5, 
C; = | 


= Ba —(5,+5,) Ca ee 
Cvs B, — ($ +5.) C) — 5,52 Cs 


C.= Bua StS Co— 55 Cara 
~ — 


4 e ; 


Cr3 — ope ~— ( StS.) 
Ch, =I 


™ 





Table 1-4, Functions of the variable coefficients f and i ° 
The functions in this table are assumed B and B, are adjustable ; 
Simultaneously. Where (-S,) and (-S ) are the arbitrarily chosen 
roots and C's are the coefficients of the reduced characteristic equation. 
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C ne = Baer — Cra (5452453) Cos (56455455 
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C= Ba- G5 4S, 45,) —C(55, #55155) 
Ge) 5: /'5,53 






functions of the 
variable coeff. 


‘f 


B54 (5&4 S'S £55) '5 ic (fs Bs 4'6)°D = a 


qa 
—SS 


| £s%5'5 + (155 + £55455) D+ (€54 %54'5)'D4°9 


Table 1-5. Functions of the variable coefficients ites : 


The functions contained in this table are assumed B 
adjustable simultaneously. Where (-S 


), (-S)) and 


> B, and 
-S are 


chosen roots and the C's are the coeetiictents of the réduced 


equation. 
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St5'5 =) + (*5'S 4 SO + en). 


f,: and ae : 

B. are 

tRe arbitrarily 
characteristic 





coeff. of the reduced functions of tne 
chara. equation. variable coeff. 


A-+h = B.~(5,+5,+5;) 


ie c,(5,45+ 5) ~( SF +$,5,45,5, ) 
5 th B, —( 5,4 +S) 


Co= By ~ Cy (5,458,455) — Cul 5,5, +55 t 5357) 
or Ba C, (5, +5,+5;) —_ C; ($5, +55; t 5,51) 
C1 = B.— C3 (5, +5, #8) — Ch (5,5, +55; $5) 


¢ @ 


a 
. ° e 


| Cyua=B,, - 65,45,45) 
C3 = | 


55'S 'D + ('s's+ 8575 + 75'5 2D = 


i. 
w~ 
M 
aA 
WY 
\ 
“} 
We 
a 
NY 
nm 
WwW 
i 
3 
Aa 
uv 
+ 
uM 
WY 
ee” 
+ 
i 
Vv) 
” 
Ce 





‘Table 1-6. Functions of the variable | 
coefficients - 
The functions contained in this table are assumed B > Foot = fe 
adjustable simultaneously, Notations are the same 4s Table 57 
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CHAPTER II 

REDUCED CHARACTERISTIC EQUATION 
2-1. Introduction: The reduced characteristic equation determines all the 
roots of the controlled characteristic equation and consequently deter- 
mines the response of the system. Usually the zeros of the closed loop 
system are known, then the total performance of the system depends upon 
the choice of the poles (roots of the characteristic equation). In 
order to establish a guide for the choice of the roots, the properties 
of the reduced characteristic equation are investigated. For the reason 
of simplicity and convenience most of the examples in this chapter take the 
coefficients themselves as variables. For the case of the adjustable para- 
meters scattered in all the coefficients, the analyses are the same. In 
fact, where the analysis of the reduced characteristic equation alone is 
concerned, it does not make any difference about what type of compensator 
is used. The analyses in this chapter also assumed that no restrictions 
are imposed in the variables, consequently no restrictions are imposed on 
the arbitrary roots. 
2-2. Coefficients as variables. 

Treating some of the coefficients of the controlled characteristic 
equation as variables provides a convenient and relatively systematic 
means for the purposes of analyses. However, the actual applications 
have some difficulties. Theoretically, the variables can be transformed 
analytically or graphically from one set of variables to another by the 
root-coefficient relation and the definition of variables. Thus the 
variables which originally are coefficients of a characteristic equation 
can be transformed to another set of variables such as the actual adjust- 
able parameters, However, the transformation may be very complicated 
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because of the inherent nature of the algebraic equations. In actual 
design, it is unnecessary, because the actually adjustable parameters 
can be treated as variables directly by the partitioning processes 
described in the last chapter. 

Mitrovic’s method treats two coefficients as variables. The applica- 
tion of Mitrovic's method to the case in which the actually adjustable 
parameters enter some of the presumed fixed coefficients requires that 
transformations or approximations must be made either graphically or 
analytically. Some of the work to this end has been denen ee 
2-3, Effectiveness of the adjustable coefficients. 

The effectiveness of stabilizing ability of adjustable coefficients 
are discussed in this section. Since only a qualitative analyses, i.e., 
the tendency of stabilization is required, the root-locus method is 


applied to the reduced characteristic equation. 


Consider a fourth order characteristic equation (2-1 the effectiveness 


) 3 2 
Sear 8 Se. NO (22%) 
of adjustable coefficients are analyzed in several cases: 
Case (1): Bo adjustable. Since only B, is adjustable, then the 


roots are constrained by B., B, and BY. Assume the reduced characteris- 


a at 
tic equation is 


ae GS es CO) (2-2) 


and the arbitrarily chosen root is “Sis then 


B, = a, 


wo 
t 


Cee ior 


Ss C= 
27 





From the above equations, solve for C 


1 C, and Co» obtain 
C, —— B; cae 5 (2-3-1) 
=e B, ae Bo 
eo, — 5S, ( B,C.) es, ce (2-3-2) 
(e353) 


Substitute (2-3) into equation (2-2), one obtains 


Bo 
ha 





S'+ (8,-S)5 +(e)5 - eo 2 O (2-4) 


Equation (2-4) is the reduced characteristic equation with the coeffici- 
1° Assume Sy has been 


chosen, the tendency of the root locus of equation (2-4) by varying By 


ents expressed in terms of the arbitrary root ~S 


is to be investigated. Equation (2-4) in root-locus form is 
Bo 
“et (S > 55) 


= Br 
s[s"+(8,-s)s+ 3 


== | (2-4) 
Consider equation (2-4) and assume S,7=0, then there is a zero in the 
right-half plane and the angle-criteria of the root-locus is 2/.M (n = 0, 


1, 2,) instead of m/7Ut (m=1, 3, 5...). The general root-locus pattern 


is shown in Fig. 2-1, in which Ry and R, are the roots of the equation 


JG) 





Fig. 2-1 Root-locus of equation (2-4) Case of 
Bo adjustable, 
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Case (2): B., adjustable. Now consider the same fourth order 
characteristic equation (2-1) in which only the coefficient of the 3rd de- 
rivative term is adjustable. By the same manipulation as Case (1), the 


coefficients of the reduced characteristic equations are obtained as 


follows: 


Cc, = B,/S, B./s; DY) 
Cy Spey Sn B/S (2-5-2) 
Gy ba, S, (2-5-3) 


Substitute (2-5) into (2-2), the reduced characteristic equation be- 


comes 


2 ! - ; Bo c - —— 
S+(5,~ ISHS )S FRCS - 55 tS) =O 8) 


i 


In root-locus form for BA as variable, (2-6) becomes 





Bo 

5? ( ie 5,5 + 5, J 

i Bo By i ae =| oa) 
ta aa a 


Assume Sy ee 0, equation (2-7), has two complex zeros in the 
right half plane and the angle-criteria of the root-locus is + 180°. 


The root locus pattern is shown in Fig. 2-2. 


Fig. 2-2. Root-locus of equation (2-7), B 
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Case (3): By and B., adjustable. 


Consider the same fourth order equation but with B. and BL adjustable. 
Since there are two coefficients adjustable, the reduced characteristic 
equation is of order 2. Follow the same derivation, the coefficients of 


the reduced characteristic equation are: 





! 
Cre ( Baa Ce, (2-8-1) 
Bo 
Co= 5S, (2-8-2) 
where -S, and -S, are arbitrarily chosen roots and S, + Ss, “0; S85 SE 0 


are assumed. The reduced characteristic equation is 


2 f Oo o 
S + sre (8,-55-22)s + =o (2-9) 


5,2 


In root-locus form for BY as variable is 


Be. C is 
Brissy | S— (514 5.) | 


[ ba 
ce + | 





—— | (2-10) 


This equation again has one zero in the right half plane and the angle- 
criteria of the root-locus is 2/i." (m= 0, 1, 2,...). The general 


root-locus pattern is shown in Fig. 2-3. 


Fig. 2-3. Root-locus of equation (2-10). By and 
B. adjustable. 
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Case (4): B, and B., adjustable. 


In this case, the reduced characteristic equation is of order 2 


and the coefficients from Table 1-2 are as follows: 





5 Sore 
pe een ee 2-1i-1 
€ aor S S2 a oe) Bo ( ) 
Be 
a 2-11-2 
Cc. 5,52 ( ) 


In root-locus form, the reduced characteristic equation becomes 


( 5, +2) _ 2st 
(5,52)* Bo (5 5,t>, —_——- 
ss ——— | (2-12) 
eee a) 


Equation (2-12) has one zero in the right half plane and the angle- 


eo 


criteria of the root-locus is 727M (n=0, 1, 2,...). The general 


root-locus pattern is shown in Fig. 2-4 when S, + Sy, =e). 


Fig. 2-4. Root-locus of equation (2-12). Bo and B 


adjustable. : 
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Case (5): By and B, adjustable = 


5 From Table 1-1. The coefficient 


of the reduced characteristic equations are as follows: 


C= 5 bea (2-13-1) 
Bo 
Co = Sz a: 
The reduced characteristic equation in root-locus form is 
— @ 
Seo 
a ee ee ES | (2-14) 
$1 S+ (B,-5,-5.)] 


The general root-locus is shown in Fig. 2-5. 


jw 


tO 


Fig. 2-5 Root-locus of equation (2-14). By 
and B. adjustable, 


Case (6): By adjustable. In the same way, the reduced characteristic 


equation in root-locus form is obtained as follows: 


Bo 


S 
- sa A =) (2-15) 
Ses + (B,-5,)S + (B,+5,/- 8551) 





es a 


The general root-locus pattern is shown in Fig. 2-6. 
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Fig. 2-6. Root-locus of equation (2-15) By 
adjustable, 


Case (7): B Bo and B, adjustable: 


1? 
In this case three roots are arbitrary, only one root is constrained. 


Assume S), is the constrained root, and define Co = Sys then 
— SHAN: Co 
C.= BA 515253 


The reduced characteristic equation is 


4 Bo - 

ee a 
If no restrictions are imposed on the adjustable coefficients, then for 
S 


any choice of S and 8, (Ss, 0, Ca OF ae 0), the system can be 


PZ 
stabilized for any value of Bo 
From the above analyses and the general pattern of the root-locus 
three conclusions can be drawn: 
(1) If the adjustable coefficients are not in the lowest order 
or not in sequence, as the cases (1), (2), (3) and (4), right half plane 


zeros are introduced in the root-locus plot. Consequently, the stabiliz- 


ing ability is less than that if no right half-plane zeros would be 
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introduced. 
(2) If the adjustable coefficients are in the lowest order 

and in sequence, as the cases (5), (6) and (7), no right half-plane 

zeros are introduced. Consequently the stabilization is more effective. 
(3) The most effective scheme is the case where all of the 

coefficients are adjustable as case (7). The most ineffective scheme 

is the case where only the highest coefficient is adjustable as case (1). 

In general, if there are "r'' unadjustable coefficients which have an 

order lower than the highest order adjustable coefficient, then the 

reduced characteristic equation in root-locus plot for BY as variable 

has "rr" zeros in the right half plane. Because of this property, it seems 

that if the adjustable coefficients of a characteristic equation are not 

the lower order and not in sequence, the stabilizing tendency of these sets 

is less than for the case of adjustable coefficients which are of lower 

and B 


order and in sequence. Consider case (5), in which B are adjust- 


1 2 


able, the reduced characteristic equation is of the second order and there 
are no zeros in the right half plane, consequently the system can always 
be stabilized for any value of K. However, for cases (1) to (4), it 
cannot always stabilize the system. In general, for any order system 

if the reduced characteristic equation is of second order and all adjust- 
able coefficients are of lower order and in sequences, then the system can 
be stabilized for any forward gain K. However, this is not the case if 
the adjustable coefficients are not of lower order and in sequence. More- 
over, the more zeros introduced in the right half plane in the process of 
derivation of the reduced characteristic equation, the less the tendency 
to stabilize the system. This implies that the more the unadjustable co- 


efficients which are of lower order than the highest order adjustable 
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coefficient, the less the effectiveness of the control parameters. This 
is the inherent property of algebraic equations. [t is because of this 
property that many approximations can be made by ignoring those control 
parameters which affect only the coefficients of higher derivatives. If 
the design of a control system is based on the characteristic equation, 
those coefficients of lower order derivatives should be considered first 
if it is possible. 

This analyses also explains the well known fact that a cascade filter 
is more effective than lower order derivative feedback, especially for 
high order systems. For a cascade filter, the filter parameters areé 
scattered in all coefficients of the characteristic equation in the nature 
of algebraic equations, while derivative feedback enters some of the co- 
efficients only. This analyses also explains the fact that tachometer 
feedback is good for low order systems, but may even hurt the performance 
of high order systems. 

2-4. Root-Region of Stability. 

The reduced characteristic equation is a function of the arbitrary 
roots. Mathematically for any choice of the arbitrary roots, there ex- 
ists solution of the equation. However, in control systems, the roots 
-of the characteristic equations are required to be in the left half plane, 
this requirement therefore confines the arbitrary rocts in a certain region 
in the left-half plane. This region here is called "region of stability" 
and may serve as a guide for the choice of the arbitrary roots. 

If the number of variables is one, this region is a line or lines 
in the S-plane (root-locus). When the number of variables is more than one, 
this region is an area or several areas on S-plane. The boundary of this 


region can be calculated by Routh's criterion subject to the reduced 
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cheracteristic equation. The pattern of 4th and 5th order equations for 


fa > ; ie ~f; and tT, - f3 as variables are shown in Fig. 2-7 


and Fig. 2-8. 
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Fig. 2-7 Root region of a 4th order characteristic equation 
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Fig. 2-8 Root region of a 5th order equation | 
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In Fig. 2-7 and Fig. 2-8, if a pair of complex roots are in the 
stable region, then all the roots are in that region, If a pair of 
complex roots are in the unstable region, then there is at least one root 
in the right half plane. If a pair of complex roots are on the stability 
boundary line, then there are at least one pair of complex roots on the 

yo axis or one root at the origin, The arbitrarily chosen roots, 
therefore, must be confined in the stable region. The calculation of 
the boundary is illustrated in the following example. 


Example 2-1: Given a fifth order equation 


6+ 16S*+ 12385 t 5208 + 13005 +2000 =O = (2-16) 


Case 1: Assume ae and fe are variables, then the controlled 


characteristic equation becomes 


5* + 16s*+ f.5° + £,S° + [3005S + 2000 == (2-17) 


From Table 1-2: 








Ce a(S, +5, ) (2-18-1) 
| - 

C,= 35, | 8. -Ceo (5, +52) | (2-18-2) 
Bo 

ey os. (2-18-3) 


where -S, and “S, are the arbitrarily chosen roots. Assume these roots 
are confined on a constant Jf line, then 


Se Sy — 2°), Spohn — WO, 


3 
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Substitute into (2-18) 


a B, ~2f Wn 





(2-19-1) 
nie 2§ Bo 
Ge 1), | 6, or (2-19-2) 
Com = (2-19-3) 
Wy 
where By = Moy By = 2000, By = 1300 from the given equation (2-17). 
The Routh's criterion for the 3rd order reduced characteristic 
equation is 
Ceo (2-20) 
Substitute (2-19) and the values of Bys By and By into (2-20), 
and manipulate, the Routh's criterion in equality becomes 
64 § 2 
——__——. sf 
2.6 FW, t+ oo), one 1) (222) 


Let the left side of the (2-12) be fs , then 


he pufeiRiitawdt 64 $_ 
J 2.65 Wr o, 


The 7 versus (,),, plot for given value of f is shown in Fig. 2-9. 


| 


+ [8.8 +8f 


oe a UW, 


Fig. 2-9 Graphical solution of equation (2-21) 
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From the inequality (2-21), if WW = Wy, =, , then all the roots 
of the reduced characteristic equation are on the left-half plane. Cu, 
and WJ, define the stability boundary of the arbitrarily chosen roots and 


are the solution of the following equation 


C4S ra 
2.6 5Wn t > 18.&8+85 (2-22) 
Solutions of (2-22) for O<~f <I[ are tabulated in Table 2-1, and 


the boundary for root region are plotted in Fig. 2-10. 


Table 2-1. Solutions of equation (2-22) 


S , WW), 
1 a.7 6.62 
0.9 3.28 7.53 
0.7 ae 10.5 
0.6 ZA05 12.9 
0.5 1.7 14.3 
0.4 1.35 17.9 
0.3 05 23.95 
0.2 0.65 36.1 
0.1 0.35 72 


Case 2: i , 1, variable 

Since the coefficients of the reduced characteristic equation are 
expressed by the arbitrarily chosen roots, the mapping contour can be 
chosen arbitrarily. In order to demonstrate this, let the two arbi- 


trarily chosen roots be on a constant O_line, then 


ee (2-23-1) 
$5. = atw?* (2-23-2) 
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Fig. 2-10 Stable root region of a 5th order equation with 
variable of example 2-1. 
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Substitute (2-23) into the coefficlents of the reduced characteristic 


equation in Table 1-2, one obtains 


C, == 16-20 (2-24-1) 
C,) = £20-(1b-2a)2a —-(aA+W*) (2-24-2) 
on 2000 /( a + u)*) (2-24-3) 


The Routh's criterion is C,C, »S Co Substitute (2-24) into the 


above inequality and manipulate, obtain 
(16-20) | 128 -(1b-2r2a -(a tw) [(X tW*) > 2000 (2-25) 


The boundary value of (J for the inequality (2-25) are solwed by the 
same procedure as Case 1 and tabulated in Table 2-2. The stable root region 
is plotted in Fig. 2-11. 


Table 2-2 Stable root region of a 5th order equation 
with f, and £, variable. 


OWN Stable region of CO 
0 0.975 < W < 11.29 “0.975 > Ws>-11.29 
2 “8.5/7 < UW < 48,57 
4 -6.62 i +6,62 
6 ~5.9 =n < +929 
Z -6.3 <ttw< 46.3 
7.5 -5.9 <= wW< +6,9 


Case 3: i ee variable. 
In this case, from Table 1-3, the coefficients of the reduced character- 


istic equation are: 


Ge St) (2-26-1) 
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Fig. 2-11. Stable root-region of a 5th order 
equation with -, and f, variable 
in example 2-1. 
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f 
C,= ere. [3.- Gy ese F Se Se OO 


i. — B./ $52 


Assume Sy and Ss, are on a constant 5 line, Then 


SS = 27 Wy , § =O 


e 


ae). 





SS. 


(2-26-2 


(2-26-3) 


(2-27 -1) 


(2272) 


Substitute (2-26) and (2-27) into the Routh criterion inequality and 


manipulate, obtain 





(16-26 )( 520 — £2 -IbW,, +29" JW, —= 2f§ « 2000 


J 


(2-28) 


This inequality is a function of GJ, UP to the 5th order, the solution 


is more easily found by a graphic method. The function of the left side 


of (2-28) for § as a parameter are plotted as a function of (,, , 


boundary values of Cue are evaluated as shown in Table 2-3. 


region of stability is plotted in Fig, 2-12. 


§ Stable region 
0.1 22 = A pe 
0.3 2.3 <= (UW 6 
On5 2.4 <= Wyr<7.1 
Oa Ze eee Oe LO 
0.9 2.6 <— W,< 8.6 
1.0 26SEC 27 33 


Table 2-3. Stable region of 5th equation with 
- 4 £ variable. 
3 
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Figure 2-12. Stable root-region of the 5th order equation with 
f > f variable in Example 2-1. 
' 3 
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In equation (2-26-2), = 15,220 was assumed. This implies 
$= 0. For S = 0, the coefficient Cy of the reduced characteristic 


equation is undefined. Rearrange equation (2-26-2) 
ent oa) = a Coro? 
For f = 0, §,+5, = 0 then 
Ba CO. —C; w* == © 
a 2 (2-287) 
Bw — BW 4 Be —_- O 


consider equation (even. the coefficients B,, B, and B are unadjustable 


4? “2 
coefficients of the given characteristic equation, therefore, the solu- 
tion of (2-28) defines four points ti), , and t jw, on the 
axis instead of a region. Those four points are convergent points for the 
stability boundary, i.e., if there are any roots on the yw axis for this 
type of control, the roots must be either + yw, , or + ;W, regardless 
of where the other roots are. Consider Fig. 2-13 which is the same as 
Fig. 2-12, but repeated for convenience. By the definition of stability 
boundary, if there are any roots on the boundary, there are at least one 
pair of roots on the jw axis or a root at the origin. Assume a point 
P (inside the stable region or the unstable region) is a root of the equa- 
tion, as P approaches the boundary in any manner whatever, equation (2-28Y 
implies there are either a pair of complex roots which approach the points 
oj, or t,W, , or a root which approaches the origin. In case one 
and case two of this example, the roots on the ya axis have an interval, 
and the location of the roots in this interval depends upon where the other 
roots are. In this case, the roots on the jw axis are independent of the 
other roots. In general, if the variable coefficients are not in sequence, 


then the roots on the |W axis have this roperty. In short, for this 
J P 9 
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Fig. 2-13 Convergence on the yw axis 
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type of control if the variable coefficients are chosen such that 
there are roots on the ce) axis, then those roots on the Jw axis are 
completely determined by the fixed coefficients of the given equation. 

In this example, substitute BY B, and By into equation (2-28)- 


one obtains 
W, = 2h | W,= 35.2 

The discussions above are for complex roots. For real roots the 
analysis can be carried out just by the same procedure. 
2-5. Dominant Root Region: No Zeros Close to the Origin. 

Since most of design is based upon @ pair of dominant complex roots, 
a procedure of evaluating the "dominant root region" is to be discussed. 
The definition of dominant roots depends upon the degree of dominance and 
other requirements. It can be defined as one wants. In case of no zeros 
close to the origin and if the pair of dominant complex roots are so de- 
fined such that the negative real part of all other roots (complex and 
real) of the original characteristic equation must be greater than that 
of the presumably dominant roots, then the Routh's criterion can be 


applied to the reduced characteristic equation by a conformal mapping 


process. 





(a) S-plane (b) p-plane 
Fig. 2-14 Conformal mapping of S = p-a 
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Consider Fig. 2-14(a), assume the pair of arbitrary complex roots 
(-S,) and (-S5) are the dominant roots as defined in the last paragraph 
and also assume all the other roots of the original characteristic equa- 
tion must have the real parts which are greater than "a" (a is positive 
real number and greater than or equal to FJ, ), then all the roots of 
the reduced characteristic equation must be to the left of the line S = -a 
in the S-plane. By the linear transformation S = p-a, the region to the 
left of the line S = -a in the S-plane is transformed to the left-half plane 
in the P-plane. Therefore the Routh’s criterion applied to the p-plane, 
implies all the roots of the reduced characteristic equation are to the 
left of the line S = -a in the S-plane. In addition, to put all other 
arbitrarily chosen roots (there may be more than two arbitrarily chosen 
roots) also to the left of the line S = -a, then all the roots of the 
original characteristic equation except a pair of complex dominant roots 
are to the left of the line S = -a. Here the line S = -a was assumed to 
be fixed if the dominant roots are fixed. In order to find the dominant 
root-region which satisfy the definition, ‘a’ must be a variable as the 
dominant roots are varying. Therefore the application of Routh's criter- 
ion to the p-plane defines a region in the S-plane if such a region exists. 

To obtain the transformation from S-plane to p-plane, the process is 
just a substitution of S = p-a into the reduced characteristic equation. 


Consider a third order reduced characteristic equation 
S 2 
ee ot CS. Cee © (2-29) 
substitute S = p-a into the above equation, one obtains 


(P-4)°+C,(p-a) + C,(p-a) +C, =O 
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manipulating, it becomes 


3 ws 3 2 >) 
Mec, 30) P + (C,-2ac, t3a)P+(C-CatGa-a) (2-30) 
== ©) 
Equation (2-30) is the transformed reduced characteristic equation in the 


p-plane. Now define: 


De, Sapa (2-31-1) 
xz& 2 

fe Oy 2qC, 13a (2-31-2) 

D&C -ca+Ga%- a’ Q-31-3) 


Where D's are the coefficients of the reduced characteristic equation in 


p-plane, then equation (2-30) becomes 


5S +DS+Ds+D=oO (2-32) 


The Routh criterion of equation (2-32) is 


DD) wD, oo 


As the coefficients D's are functions of the arbitrarily chosen roots and 
the definition of 'a', the inequality from Routh's criterion defines a 
region on the S-plane directly. 

The coefficients D's for a nth order reduced characteristic equation 
and its Routh's criterion inequality have been formulated and tabulated 
in Table 2-4. Since the reduced characteristic equation of any order 
system and any number of variables can be derived, the process of evalua- 
ting the dominant root-region for the definition described above can be 


applied to any order system and any number of variables. 
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Table 2-4. Transformed coefficients D's of the reduced characteristic 


equation "a" is the mapping contour, C's are the coefficient in the S-plane 







order of the Routh's 
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Fig. 2-15 Plot of inequality <-3/. 
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Example 2-2, Consider the same characteristic equation given in 


Example 2-1, namely - 


2 ; 4 2 ! os 
Pewee + eS fT,5 +t FSB, = O (2-34) 


where By =e Gress 


characteristic equations are: 


an 128, and B= 2000. The coefficients of the reduced 


C, = 16-27W, (2-35-1) 
te) £20 32s + (Aig ') Ww, (2-35-2) 
oo 2000 / Ww, aa 


where § and UU), define the dominant roots, since there are only two 
arbitrarily chosen roots. Assume Boe Pu, > L.@., all the other roots must 
have real parts greater than or equal to QW, , Substitute equation (2-35) 
into (2-31) the transformed coefficients of the reduced characteristic 


equation become: 


ee ie LF), (2- 36-1) 


Zz = 
D = 128-6t¢fW, +O - 1) Wy (2-36-2) 


t 


yt 2 oe S peter: 
D, = 200, — 128 fw, +48f uw, - (TF -Uw, (2- 36-3) 
By substitution of (2-36) into the Routh's criterion DD, > dD, and manipul- 


ating, one obtains: 


2 ; 2 2 Z oe 
IS3L GW, 16 (286°-idn - 45 (1 12°) W, + BA? < 2046 rel) 


Let the left side be FiQe)s choose ¢ and plot -— versus (J), as shown 
in Fig. 2-15, the dominant roots (defined by a = (iJ) region which satisfy 
the definition for F = 0.3, 0.5, 0.7 are evaluated from Fig. 2-15 as 


follows: 
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° = 0.3 wee UC), 6.3 


f 


@ =0.7 19 < WwW, < 4.8 


tt 


0.5 1.65 W,< 5 


This region is plotted in the S-plane as shown in Fig. 2-16, in which 
the stability boundary is also plotted, 

In the dominant root-region, assume the dominant roots are chosen to 
be -S,, -S, = 1.89 + ice (CO ees, CS) 
the secondary roots evaluated from the reduced characteristic equation 
are: 


Bees, = -2.8 + } 1.0, ( §# 0.94, G) = 2.82) 


3 4 
“5. = -6,6 
It shows that the secondary roots have real parts which are greater 
than that of the dominant roots. 
If the dominant roots are chosen as 
“S,, “S, = 2.5 + higeno2s ( P=0.5, W,= 5) 
the secondary roots are: 


*S,, -S, = 3.56 + J 2.76, ( 9 = 0.79, W,,= 4.5) 


¢2 


“Ss. = -3,88 

It shows again the real parts of the secondary roots are greater than 
that of the dominant roots. 

In the above two choices of the dominant roots, all the secondary 
roots also have greater damping ratio ° » However the dominant root region 
evaluated by this procedure doesn’t imply this restriction to the secondary 
roots, but this restriction can be obtained by choosing large value of "a" 
relative to f Wnyand by inspection of the shape of the region. 


In Fig. 2-16, it can be seen that the dominant root-region is inside 


the stability root-region. The secondary roots are not necessarily inside 
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SS. UU ee) Be - a OO 


the stability boundary. The max. damping cf the dominant roots can 
also be evaluated from the dominant root-region. 


Example 2-3. Given a 6th order characteristic equation 


& cc cy. ee: ee mes 
meeee> + BS t fs + f,5 +75 + oo 


= 184, By, = 10840 and BR = ox10. 


Here three coefficients are adjustable, three roots can be chosen 


where ff, , Tz are adjustable, B, 
arbitrarily. Assume the requirements of the dominant roots is such that 
all the secondary roots must have the real negative parts at least twice 
that of the dominant roots, then the mapping comtour is a = 2 pu » In 
order to meet the requirements of the secondary root, the third arbitrarily 
chosen root must be set at least to the min. value, namely S. N26tUn . 


Assume S, = 26 Wy then the three arbitrarily chosen roots are: 


ee eee 


“>, 72 a - FW, FJ oa ie -- 7 ~ - dominant 
5 -2F Wn 

From Table 1-5, the coefficients of the reduced characteristic equation 
are: 

Ge 5. —( $1577; ) (2-37-1) 

Mo 6 StS ee ce) (2- 37-2) 

2-37-3 
C, = B/ SX 5; ( ) 
Substitute S,, S, and S, into (2-37) and manipulate, obtain 
5 2=38-1 
Cs — De 4°0n ( ) 
2 2 
Gy — oa a BW, + ¢ 12° -t) Wy (238-2) 
3 
Co = Bf 2h ur (2~ 38-3) 
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Dominant root region of a sth order equation with fJ,, f, variable, 


Fig. 2-16. 
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Substitute (2-38) into the coefficients of the transformed reduced character- 


istic equation in Table 2-4, obtain 
iy =o oe [2 PU), 


D 


= 2 os 2 3 
D = B /29w, ~2§8,W, + l2f BW, - 29 (246-1) Wy 


The Routh's criterion DoD, > D, after manipulation becomes 


3 2 2 2 : 
Ww, | 8,8, ~ 8F(B+B)w, + 8, (oBs™ wh - BPl4ep~ YW | > B./2¢ 
B 


8. > Sf Bon + (40¢°- 1) Wy 


Substitute the numerical values of B 4 and Bas it becomes 


5? 


xr i. 7 2 3 hexio® 
Wy LS6XI0 -$§ (42) 10°, 4 184 (108 - )w, - §F (44§ I), | = Ff (2-39) 





For y = 0.3, 0.5, 0.7 and 0.9, the region of (J), which satisfy the in- 


equality (2-39) are evaluated and plotted as shown in Fig. 2-17. 


e = 0.3 7.5 £W, £30 
ep = 0.5 7 <5 W419,5 
a 0.7 boat) = 3 

ae Oc None 


Within the dominant root region, if the dominant roots are chosen as: 


§,, -8, = “9.5 + j 16.4; ¢ i = 0.5, W,= 19) 


eGur= -2 CLV =< 
then: S3 2 a nN 19 


The other three roots evaluated from the reduced characteristic 
equation are: 
-S = -20 + 1253 ( , = 0.97, UW, = 20.6) 
“Se = -106 
If the dominant roots are chosen as 
soy, 55 = -8 + j13.8, ( (F=0.5, WwW = 16) 


The secondary roots are 


a9 








3 
“Sys “S. = -23 + fee ( ys O87, = 26.4) 
“S¢ = -106 


The evaluation of the dominant root region is not necessarily carried 
out point by point, only those points of interest are evaluated. And 
furthermore, only the boundary points which satisfy the Routh's criterion 
are of significance. The Routh's criterion formed in this process is a 
polynomial of ),, in any case, therefore for higher order reduced 
characteristic equations there may be more than one region which satisfies 
the criteria. But only the region which is closest to the jW axis is of 
interest, because the roots are assumed to be dominant. 

2-6, Modified Dominant Root Region - One Zero Close to the Origin. 

The assumption that the response of 4 system can be characterized by 
a pair of complex dominant roots is based on the approximations that all 
other roots are either too far to the presumable dominant roots or close 
to a zero. Consider Fig. 2-18, the roots (-S,) and (-S,) are said to be 


dominant 





Fig. 2-18 General Closed Loop Pole-zero Configuration. 


60 





a re > ee aa, Oe 


because the other roots (-Se)s (-S¢) and (-5,) are far to the left 
and (-S,) is close to a zero (-2,) so that they contribute negligibly to 
the transient response, [f no zeros lie to the right of the dominant 
roots (close to the origin}, then all the secondary roots are required to 
be far to the left if the presumed dominant roots are dominant. If there 
is a zero close to the origin, then one root to the right of the domin- 
ant roots has negligible effect on the dominance if that root is close enough 
to the zero. Therefore, the dominant-root region defined in the last sec- 
tion must be modified when there are zeros close to the origin. Usually, 
the zeros of a system are known. If the known zero is far to the left of 
the dominant roots, then the dominant root region defined in the last sec- 
tion is still valid. When the zero is known to be close to the origin, 
the modification of the dominant root region is required. 

Consider Fig. 2-19(a) which is a third order system with a cascade 
compensator. Assume the compensator is a lag network. Let (-S,) and 
(-S,) be the arbitrary roots (also the dominant roots in this case), then 


(Z) and (p) are functions of S, and Sy (Equation 1-22). The general pattern 


1 
of the root locus is shown in Fig. 2-19(b). Since Sy and Sy are arbitrary 
roots, then from the tendency of the root locus, the one only root which is 
possible close to the origin is (-S,). {Tf the arbitrary roots are chosen 

within the dominant root region (point B in Fig. 2-19E) defined in the last 
section and the pole (p) and zero (Z) of the compensator are adjusted ac- 

cordingly, then (-S4) and (-S,) are to the left of the dominant roots (Fig. 
2-1Sb). Lf S, and S, 
gion (Point C in Fig. 2-19E), then the root (-S,) is forced to be in a 


are chosen on the boundary of the dominant root re- 


position as shown in Fig. 2-19C. As the arbitrary roots move toward the 
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Fig. 2-19a Cascade filter. 
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Fig. 2-19b Root locus of Fig. 2-19(a) in case Z>p 











Fig. 2-19C Dominant roots on the boundary Fig. 2-19D Dominant roots 
of the dominance for a = fW, on the boundary 
of stability 
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Fig. 2-19E Stable and dominant root-region. 
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Fig. 2-20 Dominant root region (a = fW,) and Modified Dominant 
Root Region 
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origin and so does the zero. When the arbitrary roots reach the 

boundary of stability (Point D in Fig. 2-19E), (-S,) reaches the origin 

and the zero and pole of the compensator also reach the origin as shown 

in Fig. 2-19D. From this analysis, it can be visualized that if the 

dominant roots are chosen close to the stability boundary and within 

the stable region, this implies a dipole close to the origin. Therefore, 

a region which is close to the stability boundary and between the boundaries 

of stability and dominance (which are defined in the last section cor a= 
F Wr ) is also a dominant root region. This region is given the termin- 

ology "modified dominant root region" in order to distinguish it from the 

definition of the last section. These two dominant regions are shown in 

Fig. 2-20. 

The pole-zero configuration discussed above is often the case for 
systems of high velocity error constant (K). For very high velocity 
error constant, a dipole near the origin is necessary. 

The analysis above is based on the tendency of the root-locus and 
the known zero. The same anattve te can be applied to other cases. 

In this case the dominant roots should be chosen either close to the 
stability boundary or within the dominant region (a = §Wy ). Choosing the 
dominant roots in the vicinity of the dominant boundary (a = §Wn ) 
should be avoided when a zero is close to the origin because the root is 
not close enough to the zero and it may cause excessive overshoot. If a 
zero is known to the right of the dominant roots, a root close to it is 
necessary. 

Now consider the pole-zero configuration of Fig. 2-19D which is the 
situation of stable limit. In this condition, the pole and zero of the 


compensator are at the origin, therefore the presumed dominant roots “5, 
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and “Sy are on the original uncompensated root-locus (K as variable). 
This provides another means to evaluate the stable boundary, namely the 
stable boundary is the uncompensated root-locus itself in this case (lag 
compensator). This boundary line is one of the dividing lines of the 

root relocation zone discussed in reference 7. For high order systems 

J 

there are several branches of the stable boundary line, but for the pur- 
pose of dominant region only the one which is closest to the j4) axis is 
of interest. From the root-locus it can be visualized easily that the 
dominant roots are within the region which is the closest to the jWaxis. 
From this analysis, it concludes that for lag compensator the region 
closest to the nearest branch of the uncompensated root-locus is the 
dominant root region. And the dominant boundary of the definition a = PW), 
provides a guide of how close they should be, 


Example 2-4, In Fig. 2-21, K is not to be reduced. Find the dom- 


inant and modified dominant root region: 









© 
S(s*+2541.2) 






A(S+ EZ) 
CS ey. 


me em a er er lt A NT NN mm 


Fig. 2-21 Example 2-4, 


The expressions for the dependent variables Z, p, on and C, have 


been derived in (1-22) and (1-23). Substitute S, +S, = eae, , SS, = 


CJ , B = 6, B, = 1.2, Bz = 2 into (1-22) and (1-23), one obtains 


n oO 1 2 
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ee (1-49 we - 1.2 $4.0, | b /( 2W, eeu, a) @a40-15 


p= WC, / © (2-40-2) 
Se 2 TP - 25 Wh, (2-40-3) 
a Onc, t2PWACy - L2P] /6 (2-40-4) 


From Table 2-4: 
(rec 


WRG Ca + Oo 


O 
For a = GW), 

D, = C,- 25%n (2-41-1) 

), i. SUC, i) Ws, alse) 


The Routh's criterion for D., and D , is D. => 0 and D7 0. 


1 1 


The stability boundary line (root-locus of the uncompensated system) 


is shown in Fig. 2-22a, In order to evaluate the a = (W,, dominant 


boundary line, Cand C, must be calculated. For < 


results of computation from equations (2-40) are tabulated in Table 2-5a 


0.3 and f= 0.5, the 


and Table 2-5b respectively. The dominant boundary line for a =fu,, is 


Table 2-5a. Example 2-4 Lora, = 0.3 


Wn Co C, p oN z D De 
0.4 0.65 1.777 0.01735 0.0698 0.0249 + + 
0.5 0.473 leet Z 0.0197 0.0914 0.216 + + 
0.6 0.278 1.657 0.01662 0.118 0.141 + = 
Osa 0.0539 1.584 0.0044 0.1322 0.0335 + - 
0.8 -0.085 
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Table 2-5b Example 2-4 for y= Oe 


C)) Co c p Of Z aD ee 
0.2 0.81 1.8 0.0054 0.0543 0.0995 + + 
O.3 0.615 1.7 0.00923 0.0545 0 1695 + r 
0.4 0.418 1611 0.01115 0.0685 0.1625 ~ + 
0.5 O7213 1.509 0.00886 0.079 Ort 2 x - 
0.6 0 


Shown inFig, 2-22a, If the dominant roots v = 0.5, (J= 0.55 are chosen 
which is close to the stability boundary, the closed loop pole-configura- 
tion is shown in Fig. 2-22b. For this choice of the dominant roots, the 


compensator is 
0.084( 5 +0066) 
S + co00543 


Ge 
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Fig. 2-22a Dominant root region of Example 2-4, \\ 
Shaded area is the modified dominant ° jQ) 
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Fig. 2-22b Closed loop pole-zero configuration 
of Example 2-4, 
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CHAPTER IIf 
FUNCTION OF VARIABLE PARAMETERS 
3-1 Realization Problem: 

The result of the partition process employed in this method is to 
separate the dependent variables into two parts: the coefficient of the 
reduced characteristic equation and the variable parameters (parameters 
of the compensator). The design is carried out essentially by two in- 
dependent procedures, One is to choose appropriate closed loop pole- 
zero configuration from the reduced characteristic equation and the other 
is to realize the compensator according the chosen pole-zero configuration. 


Consider a cascade compensator of the following transfer function. 


K.( 5 + aS +Q) 


: Ga) 
Seto as 


G.= 
If there are no restrictions to the type of compensator, the parameters 
ay» a b, and bo are linearly dependent in the root-coefficient relations. 
Their solutions from the partition functions are always real. For the 
function of (3-1) with all real coefficients always can be realized either 
by a passive network or an active network. However, if the poles of the 


compensator are required to be on the negative real axis, then (3-1) be- 


comes 


G = Ka Gate ans eu) 
ee ae 


(3-2) 
St P,)C Stati ) 


From (3-2) and (3-1) 


P +P. = 5, 4 PR = be C=) 





In equation (3-3), Py and P. are quadratic functions. For certain choices 
of the roots the solution of Py and Py may be complex numbers, although 

Db, and b. are always real. Therefore the choice of the arbitrary roots 
are not only necessary in the dominant root region but also necessary in 
the region which gives a physically realizable compensator. For certain 
kinds of specifications and restrictions, the compensator may be forced 
to be an active network, but for most of the systems, the designer should 


try a passive network first. The same realization problem may arise from 


feedback compensators and even pure derivative feedback. 


3-2 Linear Dependence. 

If the variable parameters are linearly dependent in the coefficients 
of the controlled characteristic equation, then the compensators are 
physically realizable passive or active for every choice of the arbitrary 
roots, because solutions of them are always real numbers. For this reason, 
the variable parameters in the coefficients of the controlled characteris- 
tic equation should be arranged (as possible) such that they are linearly 
dependent. In many cases, the linear dependent relations can be obtained 


by transformations. Consider a single section cascade compensator 


_ eS) 


a 3-4 
Ge (54 p) (3-4) 


at K. » & and PP are variable parameters then the term Kz is a 


quadratic term. However, by defining Kz =X then (3-4) becomes 


es 


For the transfer function of (3-5) alone the parameters are linearly de- 


setae 
ee ane 





Where Ke= K/ 2. (3-5) 


pendent. For the same reason a single section compensator in the feedback 
path the relations can be linearized in the same way. 
The linear dependence of the variable parameter is a sufficient condi- 


tion for the compensator to be realizable, but it is not necessary. The 
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necessary and sufficient condition can be found. For example, a single 
section cascade compensator, the roots within the single section relocation 
zones formed by the rules of reference (7) are the necessary and sufficient 
condition. 

Realization problems may even arise for pure derivative feedback. 
Consider Fig. 3-1 in which the point 'X" is assumed to be able to insert 
a summer and the signal at which can be measured. Assumed pure deriva- 
tive feedback is devised as shown in the Fig. 3-l. From rules of signal 
flow graph the determinate of the system contains a term which is the pro- 


duct of the two nontouching loop transmissions. This term contains the 


———_ eS a 
V (S$+4)(S+b) O S (see 


h, S : h,s 


Fig. 3-1. Pure Derivative Feedback 
product of h, and h,- For a certain choice of the roots, hy and hy may be 
complex numbers. For some purposes in control systems, the system may be 
expected to be overdamped and to operate in a given set of roots such as 
discontinuous operation in reference 9. For that purpose the scheme of 
Fig. 3-1 should be avoided. For pure derivative feedback to achieve a 


given set of roots certain rules and methods have been worked out. They are 


discussed in detail in Appendix I. 
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3-3 Choice of Arbitrary Roots. 

Any root may be chosen as the arbitrary root. For convenience 
the dominant roots are usually chosen as the arbitrary roots. If there 
is only one arbitrary root the damping ratio . is the convenient in- 
dependent variable. If the arbitrary roots are more than two, two of 
them are chosen as the dominant roots, others may be either fixed at a 
certain value or be defined to have a certain relation with the dominant 
roots. Suppose the original system has a zéro which is to the right of 
the expected dominant roots. Then in the case of three arbitrary roots 
one of them may be chosen to cancel that zero or be close to it. In 
general, the more the arbitrary roots, the more the design freedom one 
has. The forward gain K may be always treated as a variable parameter, 
so one additional arbitrary root may be introduced even though K is us- 
usally specified within some limits. By treating K as a variable para- 
meter, sometimes it may simplify the relations, but sometimes may make 
the problem more complicated. It depends on the specifications and the 
nature of the compensator, but the increase of design freedom is always 
Crue. 

For the case when a large number of arbitrary roots i.e., a large 
number of variable parameters is involved, a clear picture of the closed 
loop pole-zero configuration must be assumed. With that pole-zero configura- 
tion the arbitrary roots can be properly related to the dominant roots. In 
the mapping process from the contour of the arbitrary roots to the vari- 
able parameters, only one independent variable is varied, i.e., the para- 
metrical line is obtained in the parametrical space (or coefficient space 
of the reduced characteristic equation). There are infinitely many para- 


metrical lines in the parametrical space, but only those which are within 
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the specification limits are of interest. 


3-4 Coefficient as Variables: Mitrovic's Method. 

In case of coefficients as variables, this method provides a relatively 
simple relation. For lower order systems and for derivative feedback the 
application is direct. For the case of any two coefficients taken as 
variables, the function of the variable coefficients are the same as Mitrovic's 
coefficient functior which has a certain interpretation from the mapping pro- 
cess, In order to show the coefficient functions are the same as that of 


Mitrovic's method, take a 5th order equation with Bo and B, as variables 


1 
from Table 1-4, and make substitution of 51 +5>5 = 2° Wy and S485 = WW, 
from which one obtains: 
Cc, = B,— 2fWy 
y 2 
eS — Be 25 Wn ee 2 FW) = Wy 
z 2. 
= B,- 25 BW, +(C4f-t)Wy 
z u 3 
G- = oar 2° B, Wy, +(4°-1) BW, - [2¢(4¢=1) -28 | Wy, 
< ru 2 2 3 
fax W,[ 8-278, + (VGH - 280%] Bey 


2 2 3 2 eee 
$, = rF0n [B,- 29 By t (4H) By ~ 26 (HP, ] + (6-276. +f Jo 
2 8 we 429 (Ye HI) BW, + (4P-lo” 
= 24 6.0, + (1-4F JB, Sf ef qn T (3ax 
Equations (3-6) and (3-7) are the same expressions as obtained from Mitrovic's 


meeeod: By the same manipulation, functions from Table 1-1 to Table 1-3 
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are the same expressions as obtained from the modified Mitric’s mecno die 
Mitrovic's original coefficient function is derived by mapping a constant 

‘ contour through the characteristic equation. By the same mapping 
process, other contours also can be used. In this analytical method the 
functions are expressed in terms of complex variables, the functions for 
any contour can be obtained readily. Because of this extension of Mitrovic$ 
Method to the contour other than constant , line, some of the features 


are discussed briefly in the next three sections. 


3-5 Cauchy's Theorem Applied to the Function cf Variable Coefficients. 

In the case of two variable coefficients of a characteristic equation, 
the interpretation of the parametrical line of f VS ie is discussed 
in detail in reference (2) and reference (4). Here only the cases of 
variable coefficients which are not in sequence are discussed. Take the 
case of a and fe: as variable coefficients. Assume the original character- 


istic equation is 
n 3 2 
Se—ooete- -- + 6,5 + B.S + BS + Bz (3-8) 
and assume the controlled characteristic equation is 


3 2 
Ge ---- i fs + Bs + F548, (3-9) 
If f, and , are so chosen such that £ = 0, then 
3 
F = (8,-t)s?+(8-f)S (3-10) 
Assume the arbitrary roots are on a constant f line in the left half of 


the complex plane, then the angle of the 5? term of (3-10) is YC and 


the angle of the S term is /P if B. > ice and By an respectively. 
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rf B. <= fe and By < ie , then the angle of the corresponding term 
is /-36 and /-© respectively as shown in Fig. 3-2(b). The functions 
of +, and f, have negative values as W<<{ from Table 1-3. Assume all 
the roots of equation (3-9) are to the left of the F line, then the mapping 
of W, from 0O-co line, must give counterclockwise loops. At the 4 points 
on F(5) curve as shown by (1), (2) (3) and (4), f, and i; must have 
the following values: 

at point (1) aie aro: een 

at point (2) i = be ee (o> 

at point (3) i, = B3 oe 

at point (4) i red B, +, = 8B, 

By tracing the curve + vs i by (3-11), the it vs f, curve must have 
a counterclockwise loop and enclose the point M. In general, the curve te 
vs F in which | 2+ 1, then a counterclockwise loop of the curve ie vs T; 
implies all the roots are to the left of the F line, i* point M is 
inside the loop. In addition to the conclusions from the reference (2) and 
reference (4), the interpretations of vs i curve are summarized as 
follows: 

(l) If j=i+1(j2+0, 1, 2,...), clockwise loop of a vs i 
curve implies all the roots are to the left of the § line providing point 
M is inside the loop. 

(2) If j> i+1 (j= 0, 1, 2,...), counterclockwise loop of f. 
vs tT; curve implies all the roots are to the left of the re line pro- 
viding point M is inside the loop. 

The above rules can also be applied when the mapping contour is a 


constant vertical line ( O constant). 


76 





(4) S—plane 


(6) Fis) 
Plane 


(C) fF, vers f, 
plane 





Fig. 3-2 Mapping of constant 9§ contour to coefficient 
| plane 
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3-6 Evaluation of roots from the i vs f; plot. 

The 1. vsj;plot is a mapping process from the complex plane. The 
mapping functions of this process are the functions of the variable co- 
efficients, namely ie = (3) In this mapping process regions map 
into region, curves map into curves. An unbounded region on the S-plane 
maps an unbounded region unto the coefficient plane, while bounded regions 
map unto bounded regions. Therefore, the stable root regions discussed in 
the last Chapter map a stable region in the coefficient plane. The map- 
ping process from S-plane to coefficient plane is unique, but the inverse 
(from coefficient plane to S-plane) mapping is not necessarily unique, it 
may have multi-values. By using this multi-value property of inverse map- 
ping, more than one pair of complex roots can be evaluated from the coef- 
ficient plot. Figure 3-3 is a 4th order equation with fe and ie as vari- 
ables. Fig. 3-3(a) shows the stable root-region and Fig. 3-3(b) shows the 

af vs. i plot. It can be seen from Fig. 3-3 that the unbounded stable 
root region on S-plane map an unbounded region unto f ~ a plane. 

The s = 0.5 lines on S-plane map a clockwise loop unto the plane. 
Every point on S-plane map a point on f,»f, plane. However, the inverse 
mapping is double-valued, such as points (2) and (3) on oa plane 
map two points (2)' and two points (3)" on S-plane. Therefore, if suf- 
ficient curves have been drawn on f.~ f; plane, by using different con- 
tours (such as constant § and constant @ lines), two pairs of complex- 
roots can be evaluated from the eet plane directly. For the case of a 
4th order equation, all the roots can be seen on the tT; ~T; plane for 
every choice of f; and ae . Real roots are not shown on the coeffici- 
ent plot, but they can be evaluated easily from the reduced characteristic 


equation. 
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Fig. 3-3 Mapping of 4th Order Equation 





Since for every region on the S-plane, there is a corresponding 
region on the coefficient plane, the required sya root-region can be 
mapped on the coefficient plane if such a regia exists. For example in 
Fig. 3-3, if the roots are required such peat all roots have Oe 
and time constant less than 1, then tg root region of the equation has 
iy plot which is the intersection re- 


a corresponding region on the ~ 


a” Af , ° ° 
gion of the clockwise loops 944 = land § = 0.5. This region is shown 


by the shaded area. For. cher requirements, the region on the coefficient- 
plane which satisfies the requirement can be located by the same procedure. 
Figure 3-4 fd Figure 3-5 are numerical examples with variables a 
and a for ary/order and 5th order equations respectively. In Figure 3-4 
only one ofMstant CX (= 1) contour is drawn and in Figure 3-5, only one 


Pestagl 3 ( = 0.5) is drawn. From those figures, it can be seen that the 


constant and constant © contour constitute curvilinear coordinates 


on the mao plane. From the point of view of curvilinear coordinates, 


two pairs of complex-roots can be read for every point on the ir t- plane. 
3-7 Stability boundary on the coefficient-plane when the variable coefficients 
are not in Sequence: 

As discussed in Section Three of Chapter Two, when the variable 
coefficients are not in sequence, the functions of the variable coefficients 
are undefined on the yw axis. Consequently, the stability of the bound- 
ary line on the coefficient-plane cannot be obtained by mapping the “Waxis 
of the S-plane. However, the boundary on the S-plane is well defined as 
has been shown in Chapter Two, therefore, the boundary on the coefficient- 
plane can be obtained by mapping those boundary lines from the S-plane to 
the coefficient plane. This process can be carried out with point by point 
calculations or by approximating the boundary on the S-plane by a function. 
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The boundary in Fig. 2-12 has been calculated point by point and plotted 
in Fig. 3-6. In Fig. 3-6, a constant ° = 0.1, line is also mapped onto 
the ic ie plane. The i = 0.1 line shows a counterclockwise loop which 


implies that there exists a region in which all roots have J >0.1. 


81 





= Constant Yv\ 
constant 


“O1ila ilne 
Dash line 


proc of &@ 40M OLTgder equatllon. 


I-4( a) 


rig. 


PERERA 


ch 


EEL 


2 * ak ; = 
' | 
. 7 
ace =i =| |e - 
{ 
S no + ; — sais P 

t i } 

{ me 4 ba ‘ = = 





er 
Soe 


2 
he 
meer 
fe 
Pol ste boa 
ctEP 


{ i ni 
; ' 4 t om ¢ 
we Fe ne we 


° e 
— ial onde § ames wy — 2 ee 
4 


om 
~ owed « —_ owes _— 
' ‘ 


So terre 


=| 





Fig. 3-4(b) 


ve 


Root- region 


83 


LAK 








ty “ #3 pot. 


olaDLILly ooundary on 


y-9 


Pilg. 












CECE EE reece Cee ere eee eee eee eee eee eee EEE pce EE EEE eee ESHER EE PEE EEE 
Pee ee eo BECERRA EERE EEE atte CEE 
SES SESE SSE ESSEC eS Ee Et SESH EEE 
BEECHES eee eee eee eee PE eo PCH EEE HEE Pee ones Ht cuueeueane 

BERR EEEEEEEE EEE ESC EEE E EEC CEEECEREECCEES EEE EEEE EEE EEEEEECEEEEE EEE CEEEEEEEEECEEEEE C 

HERE eee 
RS Ee re eee epee ee Ee eo eco eee ere eee eee eee eee ee eee eee eee 

fee intet cee fees Se eae eee ee ee ee eee eee eee eee eee ee eee eae irene ete 

Se ES Se Se aS ea EEE Sate apatte 

23 Fe ee SEES ectistg EEE eae ee SSSR rE EEE ef ena 
SECC SS Ee ee eee PEE EE oe tt PEE SEC pt To ro feet Pooh neces oH Ree 

SEE CE EEE reer ee eee eee eee eee ee eee eee eee eee eee eee cee ee cere ee ee eee eee ree eee eee 
sag | SS a eg atanstogssp peg doers capresso ade deli PST sat Sly cont | Popo SSIS Pe SS 
BEER EEE EERE CERES EEE EERO REESE EERE EERE Eee EECEEE SE CEES eee ee eee eee ee eee eee Ce 

SEEPS Ecrece rere ere ep ere eee eee eee eee re rere eee eee Ee eee eee 
EEE EEE ECE ECE EEE ESC EEE CECE EE CESSES ee 

HE Fer rere eee ee ee eee ee eee ee eee eee eee ree eee eee eee eee ere eee eee 
EE am EEE EE EEE EEE ES EEE CCE EER CEE EE ee eee he BE eee 
pestis p satay, faite Heat Ter Need fiat sete ect io ab fa er HEPC TS Eth Sec Sesto tata hee 
ECE CSCC OSS EE ECE CRESS ECE EERE CESSES ECE CEEE EERE CE ECE ECE EEE EEE EERE Eee 

CEE Cre re ree re ree eee ee ore eee eee eee eee eee eee eee ee eee eee ee eee eee eee Hee 
ECE CEEEE EEE EE EEE EEE EEE CECE Es PERSE ee ee Serer Pee CEE OPE eer Eee Lee ee BEECHER EEE EEE 
BEE EEE RE ett et BE SSN H SH BECHER CEE EEE EEE EEE EEE EEE EEE EEE EEE EEE CEE 

peer CE CEE Cee EEE eS Cee eRe eee eee Cee eee ee ee ee eee eee eee eee ee eee eee 
SEE Cer Eee re eee eee eee eee eee eee eee ee eee eee eee eee 
Sues spots snnzs tones innazsnant sansa tesa, suazs anes arazeveres stats fatal sarsvanvee ated -oucseresseates tev szevend recererezefctst eves zorssntet tezezoraessctataestaestressrerazesszeresoeezetens! (at080f08 
ESSE Freee rere ree eee eee eee eee eee cee ree ee eee eee 
ESreEeeEe ESCHER EEEE NRE Cee eee eee ee Eee eee eee eee EEE SEE REE 
EEE ee He See Paap attest seh cetera 
pp tf ed es + fama eta ro neamaae po i ss eee re Sop rat PE te A Eeaeca eed oo Coo ae pee ieee et ++ = Hi + scat HH pe ft LL 
Sane See HH HH SHeteaat bails HE oe oH seseas pH +H HH oH srezais ae Seseedie ivsat etantt Po et He oH CECE HH Hoe 
He au Beeue ere ee Saree i soe aS eH 

Ht PP ts 7 ppt bt tt try ze bt | ett Het Py fy Fy ts pe ed ate | HH PLOY See eee PO see ae } ate 

er eee eggs Se seree ge gee eee eee Seo eeeheR HE aren sities loo eae Pepe 

ee CEE Ea ripe pS ECEEEEE EERE SEER EEE ECC ee 

roy oor Pore aE Reon Pee Pe eee Eee eee cee eee 

Eee eC See Ce rere rere ree eee ee cee ee ee Eee eee Eee eee eget ret Pee eee 

f ae Gauss J Pasta taa pte peta g Aa Boa Yi ph Bs fs | rt [fan rae ea ef] Ya eae al ae a Ae |e) CO a le ee Ts 

Bll ed net te etc ent tne te tar at ee lane EGE ee eee ee eae ee eee ee ec 
soit At SCE Reali SH hl Se FERRE et RECHT Ee HE et ee 

ab el ceed (ea ete (ue coed [eater (asfseeé ete Sela Sie els eset SEE SEES SSS HSHES UTE ceo feseaoepeeeeee ETT 
Ee Pe ee eee eee eco Pe eee eee eee ee eee eee ee ee 





PART II DESIGN 
CHAPTER IV DESIGN PROCEDURE 
4-1 Specifications. 

The transient response of a control system is determined completely 
if the closed-loop pole-zero (root-zero) configuration is known, The 
contributions to the transient response by each pole and zero can be 
calculated precisely. If a pair of complex roots is dominant, then 
approximations can be made. Based on these dominant roots, the effect of 
all the zeros and all the secondary roots to the transient response can 
be visualized from the complete root-zero pontigtear tous oe With this 
relation the specifications from the desired response can be conveyed to 
root-zero configuration. As in this method, the design is based on root- 
zero configuration, all specifications must be interpreted in a set of 
appropriate configurations of roots and zeros, Frequency response is 
readily interpreted if dominant roots are assumed. The steady state 
specifications (stiffness and error constant) are contained implicitly in 
the equations and formulas to work with, then they are satisfied automati- 
cally if the final results satisfy all other specifications. However, to 
meet all the specifications may not be possible until some ad fasemente are 
made. For example, to design a system of high velocity constant (K.) a 
dipole near the origin may be required. If the original dominant roots 
are defined so that all secondary roots are far away to the left of the 
dominant roots, it may result that no root-zero configurations for this 
definition of dominant roots are able to meet the high K requirement. 

In this case, the definitions of dominant roots may necessarily be adjust- 


ed so that a dipole to the right of the dominant roots is allowed, 
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In other words, a term of long tail of relatively small magnitude in 
the transient response is allowed. Usually, this long tail term in 
transient response is avoided. But in order to achieve high velocity 
error constant by a relatively simple network, this adjustment may be 


acceptable in consideration of all the factors in design. 


4-2 Design Procedure -- General. 
In the design of control systems, many methods and tools are avail- 

able such as the root-locus method, frequency domain method, pole zero 
configurations, etc. Problems frequently rise about what procedures and 
patterns are more logical. In pure synthesis, the designer is able to 
take the specifications and proceeds logically in a step by step manner 
to a system which will meet these specifications. In these senses many of 
the design techniques are not true synthesis procedures. On the other hand, 
the design procedure proposed by eaaiteniints which takes the form analogue 
to modern filter theory is a synthesis technique. In Guillemin’s procedures 
the closed loop transfer function is first determined fromthe specifications 
and the compensation is determined subsequently. The design procedure by 
the analytic method developed in this paper is more or less close to Guillemin's 
procedure. Briefly, the procedures are as follows: 

(1) The regions of closed loop pole-zero configurations are deter- 
mined from the specifications. 

(2) The type of compensator (cascade, feedback, pure derivative, 
or mixed network) is chosen from the consideration of specifications and 
other restrictions. 

(3) Analytic expressions of all the roots and the parameters of 
compensation in terms of the arbitrary roots are obtained. 


(4) Within the pole-zero region defined in (1), choose a set of 
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appropriate pole-zero configurations. Then the parameters of the compen- 
sator are calculated from the expressions in (3). 

Procedure (1) is an interpretation or transition from the specifica- 
tions to the closed loop pole-zero configurations. Procedure (2) depends 
upon the specifications, the nature of the plant and also on engineering 
judgment and experience. These two procedures are general to all methods, 
if the design is based on transient response or time domain.. They are the 
prerequisites, Procedures (3) and (4) are concerned with the method it- 
self. For this analytic method, they are explained in detail in the 


following sections. 


4-3 Design Procedure -- Detail. 
When the type of compensation has been chosen, the design procedure 
for this methid is as follows: 
(1) Formulation of the Characteristic Equation. 

The characteristic equation with the compensator is formu- 
lated in the way such that the coefficients of this equation are functions 
of the free parameters of the compensator. In other words, all parameters 
are independent from one another. 

(2) Formulation of the Root-Coefficient Relations. 

Set the coefficients in (1) equal to the corresponding root 
coefficient. If the order of the equation is "n" there are "n" relations. 
If there are "r'' free parameters in the compensator, then "r' roots are 
arbitrary and the reduced characteristic equation is of the order (n-r) 
and has (n-r) coefficients. The root-coefficients relations are formuiat- 
ed in such a way that the coefficients are functions of the arbitrary roots 


and the coefficients of the reduced characteristic equation. 
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(3) Formulation of the Functions of the Dependent Variables. 
In the "n'' root coefficients relations formulated in Step 2, 


"n'' dependent variables and r independent variables. The ''n"! 


there are 
dependent variables are the "r" parameters of the compensator and the (n-r) 
coefficients of the reduced characteristic equation. The independent vari- 
ables are the "r'' arbitrary chosen roots. Since there are "n" equations and 
"n" dependent variables, the dependent variables can be solved and expressed 
in functions of the arbitrary roots. 

(4) Determination of Stability and Dominance. 

From the reduced characteristic equation, determine the 
stability root-region and the dominant root-region. The dominant root 
region depends upon the degree and the definition of the dominance which 
is defined by the designer from the consideration of specifications. 

(5) Evaluation of the Compensator Parameters. 

If the desired dominant roots are within the dominant root- 
region which has been determined in step (4), then choose the dominant 
roots (usually independent variables) and substitute them into the formula 
of step (3), then the parameters of the compensator are evaluated. 

(6) Realization of the Compensator Network. 

If the parameters of the compensator evaluated in the last 
step are within the specified range and can be realized, then the design has 
been completed. 

(7) Design Adjustment. 

If the value of the compensator parameters are not in the 
specified range or the desired dominant roots are not in the dominant-root 
region, then it indicates either this type of compensation cannot meet the 


requirements or the definition of the dominant roots must be adjusted. When 
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other types of compensation have been chosen or the degree of dominance 
has been adjusted, the design procedure is just a repeat of the above 
steps. 

The procedures described above are general. The first three steps 
are the partitioning process as discussed in Chapter I. For some specific 
cases, formula have been derived ae shown in the subsequent sections, then 
the design is just a computation and is started from Step (4). The compu- 
tation is best carried out by tabulating all previous data and by pro- 
gressing computation. Once those data are tabulated, all information 
can be found from the table and if some changes of the type of the com- 


pensation are made, most of the data can be utilized. 


4-4 Free Parameters of the Compensator. 
In step (2) of the last section, a terminology "Free Parameters" was 
used. It may need further explanation. Consider a cascade single section 


R-C compensator. Let the transfer function be 
See K (Stz) 
GG Sie 


The apparent parameters are Kas Z and p . Assume the system is type I 
and the velocity error constant KY is specified. If the forward gain K 
of the uncompensated system was determined by the Ke requirement and in 
the partitioning processes, K was treated as a given numerical constant 
then K. is not arbitrary in order to meet the requirement of KS In other 
words K, must satisfy the following relation: 

Kn = PA 
Thus the three apparent compensator parameters Kos p » and Z are related 
by the requirement of KS Only two of them are independent variables, the 
third one is dependent. The independent variables of the compensator in 


this sense are given the terminology "free parameters", 
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For other cases, the same treatment can be applied to the compen- 
sator parameters. In the above example, if the forward gain K was also 
treated as an adjustable parameter, then all the three apparent para- 
meters are independent because KK is also independent. However, the re- 
sult of the partitioning processes are different for the two cases. But 
the design can be carried out by either case. In the former case (K is 
constant) there are two arbitrary roots, while in the latter case (K is ad- 
justable) there are three arbitrary roots. The expressions of the para- 
meters are different for different cases, but the design can be carried 


out in both cases. 
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CHAPTER V DESIGN OF PURE DERIVATIVE FEEDBACK 


5-1 Tachometer Feedback - Second Order. 
For compensating lightly damped instrument servos, tachometer feed- 


back is commonly used. For second order cases, as all information can 


be determined precisely, only the formulas are presented. 





Fig. 5-1 Second Order System with Tachometer 
Feedback 


The block diagram of a Type 1 second order system is shown in Fig. 


5-1. The separation of the forward gain by K and K, is for the purpose 


of proper adjustment of h. Assume the specifications are given Ky % and 


5 . The design procedures é2re as follows: 


(1) Formulation of the Characteristic Equation 


The Characteristic Equation is 
5S + (at Kh)S + KiKi. =0 (5-1) 
(2) Formulation of the Root-Coefficient Relation. 
Since there is only one parameter of compensation, there is only one 
independent variable. The independent variable is not necessarily the 
In the case 


arbitrary root itself, it can be a function of the roots. 


of one arbitrary root and that root has to be complex, the independent 


variable can be chosen as either iS or WW, of the pair of complex 


conjugate roots. For second order case, the only choice is § 3 
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because (,), is fixed by the gain. By doing this the root-coefficient 
relations are: 
Q+k,h = 2fp (5-2-1) 
KK, = Wh Sone) 
(3) Formulation of the Functions of the Dependent Variables. 
Since the root of the reduced characteristic equation is the domin- 
ant root itself, there are no stability and dominant region problems. 
Here the dependent variables are Wp and h( S is chosen as independent 
variable). Their expressions as functions of the independent variable are 
as follows: 
From (5-2-2): 
CAD ae (5-3-1) 
ore =e) e290, =o = 2f IK, =O. 


»¢ IRR -& (5-3-2) 
Ra 


For the given value of Ss , h is computed. 





| 


5-2 Third Order: Tachometer Feedback. 
The block diagram of a third order type one system with tachometer 


feedback all around the plant is shown in Fig. 5-2. 


Ke 
y | SGraGw 


Fig. 5-2 Third Order System with Tachometer 
Feedback . 
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The characteristic equation of the uncompensated system is 


5 +B,S+BS +B, =o (5-3) 


where Bo =a+b, By = ab, Bo = te 
Assume the forward gain K and the damping ratio f based on the 
dominant roots are given. The design procedures are: 


(1) Formulation of the characteristic equation: It is easily 


formulated as follows: 
oe BS ti SIGS. i =O (5-4) 
(2) Formulation of the root-coefficient relation: Since 
there is only one control parameter h, only one root is arbitrary. 
This root can be chosen either the real root or § of the complex roots. 
Assume the real root is chosen and designated as “3, (S5 is positive on 
the negative real axis). Then S, is the independent variable. The re- 


1 


duced characteristic equation is of order 2 and therefore has the form 
Stes +e. — Oo 

Since the roots of this equation are actually the dominant roots, it is 

expressed by P and WJ), for convenience. Here f and (J, are not in- 

dependent variables but dependent variables. By the choice of the 


above definition, the root-coefficients relations are: 


Beet +S, (5-5-1) 

2 
Brn Keee= on, + oP Ons So 
B= W,S, (5-5-3) 


(3) Determination of the functions of the dependent variables. 
Here the three dependent variables are g§ , W, and h. In the 


actual computation, only one dependent variable must be found explicitly 
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as a function of the independent variable, the others can be computed 
by using the implicitly functions of formulas (5-5). But for the pur- 


pose of illustration, they are solved and shown in equations (5-6) 


§ = (8-5) 15/2 IB 6-6-1 
ges (5-6-2) 


Kh = B/s, + BS,-S - B, (5-6-3) 


Note the three variables [ , W,and h are expressed as functions of the 
independent variable. 
(4) Determination of stability and dominant root region: 

In this one control parameter case, the stability and dominant root 
regions are lines. These lines are the same as the root locus, by treat- 
ing h as variable and K as a fixed number. The only difference between 
them is the independent variable. In the root-locus method, the indepén- 


dent variable is h while in this method it is the arbitrary root § A 


i’ 
numerical example is shown in Fig. 5-4. The degree of dominance can be 
visualized easily from those plots. 

For this particular scheme, it indicates a max. damping. This max. 
damping can be found from equation (5-6-1) analytically because , is ex- 


pressed explicitly as a function of S$ Differentiate S with respect to 


L 
S, and set it to zero, one obtains: 


1 
i = B, /3 Cx) 


Substitute into equation (5-6-1), obtain 


0. = B/218, (5-8) 


Substitute this value of S, into § and WW) 


1 the real part of the 


n > 


complex root becomes 
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Pwr —— B /3 —— ohn 


(5-9) 
Equation (4-7) shows that if the tachometer loop gain Kh is so 
B 
adjusted such that the real root is ( - =e ) then the complex 


roots have max. damping for the given value of Bf = K) and this max. 
damping is shown in equation (5-8). Equation (5-9) shows that for max. 
damping the real part of the complex roots is equa! to the real root. 


Fig. 5-3 shows this root configuration. 


-Sim 


Fig. 5-3 Root Configuration for max. § of 3rd Order 
System with Tachometer Feedback. 


For this root configuration, the tachometer loop gain and WW, are 


obtained from equations (5-6-2) and (5-6-3) as follows: 








pie. 3Bo feo 2 sk 2 
Kn ~~ ahs g B ey (5-10) 
a 
eee Ba 5-11 
on. | = ( ) 


Equation (5-8) gives a quick look of the max. possible damping for given 


K or the max. K for a given y in a third order system with tachometer 
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30 
Wh, 
02 ; 20 
a Cor 
0.| 1O 
° 2 4 © i) 10 
> 
Fig. 5-4(a) 





[2 


Plot of Equations (5-6-1) and (5-6-2) for numerical values 
given in example’ 5-1. 


WW 
= 
a oO 
a\ 8 
= wv 
x\ ®| 5 = + 
—\ TV 4k=10¢ | J 
20 
Fig. 5-4(b) Root locus of the characteristic 
Equation 
57+ 15S + (bot Oh)S +750 =0 ee 
sie 
=HS 
hk, 1S 
SF 


#k,=0 


-Cr. -20 ~{O 
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feedback. 
Example 5-1: For the block diagram shown in Fig. 5-2, if b = 10, 


aes >, K Ky = 750; the uncompensated system has a pair of complex roots 


1 


on jw axis. For this pole-zero configuration: Bo SLD: By = 50, BY = 720. 


From equation (5-8), the max. possible damping 
3 
27x 71s0 
From (5-7) Sam = 


Assume this value of § is taken, then from (5-6-2) and (5-6-3) 


For other root configurations, and (i), as a function of S, are plotted 


i 
in Fig. (5-4)(a) and the root-locus for h as variable is shown in Fig. (5-4) 
(b). If the dominant root is defined so that the secondary root has greater 


real part than that of the dominant roots, then the portion of the root 


locus when Sy —= 5 is the dominant root region. 


5-3 Third Order: Tachometer Feedback to Interior Node. 
Consider Fig. 5-5 in which point "X" is assumed to be available to 


insert a summer. The analytic design procedures are; 


"xX fe 






Fig. 5-5 Tachometer Feedback: Minor Loop 
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(1) The characteristic equation is: 


S°+ (8,4 Kh)s?+ (B+ bKh)S +B, = (5-12) 


where Bo = atb, B, = ab, BY = KK, 
(2) Root-coefficient relationship. For the same definition 
of dependent and independent variables as described in the last section, 


the relations are 


Bt Kh = 29W,,+ 5S, (5-13-1) 
B, + bk,h = Wry + 2fW,S, Cae 
R= oe So 0O-13-3) 

oO 


(3) Function of dependent variables ( f » Ws Kh De 


From equation (5-13-1) and (5-13-2), eliminate Koh 


§ = (B,0~B, tw -b5,) / 2Wn (b= Si) 
Substitute (UJ, from (5-13-3) into the above expression obtain 
eae Bb — B, + $e - 65, 


y= AS ynlne.s, 


The other two variables can be obtained in the same way but are not 


(5-14) 


necessary. 
(4) Stability and dominant root region: For the numerical 
values given in Example 5-1, the dominant roots \ and Wr are plotted 


in Fig. 5-6 as functions of S, from equation (5-14). The degree of dom- 


1 


inance can be seen easily. 


(5) Take § = 0.5 which corresponding to S, = 20. From (5-13- 


1 
3) 
« LO n= 6.14 
From either (5-13-1) or (5-13-2), obtain 
K,h = 11.14 


oe, 





ae 


| 


Pigeeo-6 


and W&, as a function of the real 
rowt 


of the systems shown in Fig. 5-5. 


O} 


O2 
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5-4 First and Second Derivative Feedback, 

In case of first and second derivative feedback around the plant 
for a given value of forward gain, only the coefficients of the first 
and second derivative term are variable. The design may be carried out 
using two variables. But for more freedom it is better to take the con- 
Stant term also as a variable then three arbitrary roots are available. 
The coefficients of the reduced characteristic equation and functions of 
the variable coefficients have been derived and tabulated in Table 1-6. 


Example 5-2. Given the original characteristic equation, 


67+ 1os?+3is° +308 +K =O (5-15) 


The forward gain K must be greater than 1000. The gain at the stability 
limit for the uncompensated system is 84, first and second derivative feed- 
back are chosen, When three variable coefficients are used, equation (5- 
15) becomes 
a S 2 a 
owes Ff, Seatel, S + f, =O (5-16) 

where f, = K 

From Table 1-6, the coefficient of the reduced characteristic equa- 


tion and the functions of the variable coefficients are: 


Cos Tom) (5-17-1) 
f. = Cy 5,525 (5-17-2) 
fe COUPER, HISES, tuS, S, ) (5-17-3) 
i, ee Sale tS, einen wees s,s, ) (627-4) 


where (-S)), (-S,) and (-S.,) are the arbitrary roots and Co = Sy: Here 
the order of the equation is four, all the roots may be complex. There 
it- 


is only one constrained root so it is convenient to replace Co by Sy 


self. The order of the reduced characteristic equation is one; equation 
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(5-17-1) alone determines all the roots. Assume (-S,) and (-S,) are the 


dominant roots, then (5-17-1) becomes 
(Oe pets TS, (5-18-1) 


f, = Wh 5,5 (518-2) 
There are no zeros in this system, all secondary roots (-S, and “S,) are 
required to be to the left of the dominant roots. Assume (-S4) and (-S,) 
are real roots, then for the definition of dominance, the max. K occurs 
when S, er Ss, §Wry . From (5-18-1) for this value of the secondary 
roots, the conditions are: 


so 2.5 (5-19-1) 


2 


2 
Max. K = 6 Wn (5-19-2) 
For § = 0.5, the max. K = 156. This is far below the required K, then 


J 
the secondary roots must be complex. Assume f and W denote the second 


roots, then equations (5-18) become 
Gud + ew, = 5 (5-20-1) 


K= Wwe" (5-20-2) 

Notice in equation (5-20-1) there are three arbitrary variables and this 
equation alone determines all the roots. Assume the dominance is defined 
as gw, > GW) then from (5-20-1), the dominant root region is 

Pw < 2.5, which is shown in Fig. 5-7a. Within this region chocse 

‘ = 0.5, ee = 2, then from (5-20-1) the condition for secondary 

roots is Qu), = 4 

Substitute into (5-20-2), obtain K — 64/6? 


For K = 1000, obtain ‘= 0.253. For this choice of the dominant roots 
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f and f, are readily calculated from (5-17-3) and (5-17-4) as follows 
f, = 532 
f, = 270 


The root configuration is shown in Fig. 5-7b. For another choice of the 
dominant roots, the computation can be carried out in just the same way. 
Example 5-3: Given the original characteristic equation 
674 16S + 1285 + Gos t 13008 + K =O 
The forward gain K must be greater than 2000. For this value of K the 
system is badly damped, first and second derivative feedback are chosen, 


By using 


a 
an 
DC: 
1 
Ww 
Ww 


-_ 
€ 





ee ee 


Fig. 5-7a Fig. 5-7b 
Dominant Root Region of Root -configuration of Example 5-2 
Example 5-2 for § = 0.5, WJ, =x 2 
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three variable coefficients, the working formulas from Table 1-6 are as 


follows: 


ae ton 3) (5-212 
eG, (Gi Ss £5) — CSS, F534 SyS)) (5-21-2) 
ieee (5~21-3) 
re oS, 5) tC, 5,5, S, nea 


en . | (5-21-5) 
T = Cost. +5.) + C(55,4 $5,455) 4 5.5.5; 
Assume (-S,) and (-S,) are the dominant roots. The order of the equation is 
five, one root must be real. Let this real root be (-S,). Here three roots 


are arbitrary. S. may assume any value. For dominance, assume 


3 
5, = 2F Wn (5-22) 
Substitute (5-22), B, = 16 and B, = 128 into equation (5-21), one obtains: 
es Oy (522321) 
Co == 128 -49Wnl, - (14 49), (5-23-22) 
eee Ste (5-23-3) 
= CW, (it 487 + ce, (5-23-4) 
£ = 4GWCo + (IF SU WRC, +29Wr (5-23-5) 


Notice (5-23-1) and (5-23-2) determine all the roots, (5-23-3) determines 
the forward given K while (5-23-4) and (5-23-5) determine the derivative 
gain. Since there are no zeros in the system, the secondary roots are 
required to be to the Peft of the dominant roots. Let a = or » the 


transformed coefficients from Table 2-4 are as follows: 
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eee or (5-24-1) 
D, = C,-Fw,c, + Pwr (5-24-2) 


The stability boundary lines are determined by C70 and oe - : 
The dominant boundary lines are determined by D>o and D,>O 
They are shown in Fig. 5-8a. For ° = 0.5, the results of computation of 


a, Cc 


B are tabulated in Table 5-l. Examine Table 5-1, 


1? fo De and D 


when WW, is between 4 and 5 the gain K is greater than the specified 


1 


value 2000. This range of (J, is within the dominant region since both 
D, and D, are positive. The root-configuration for i) = 4, ‘i meee 


is shown in Fig. 5-8b. For this choice of the dominant roots, the vari- 


able coefficients are computed readily from (5-23-3), (5-23-4) and (5-23-5). 


n= ee O20 

oO 

f, = 1534 (original 1300) 

f, = 576 (original 520) 
For other values of S., and other choices of the dominant root the 

computation is the same. 

On Cc, Co + D, Do 
I 12 V2 57 + + 
3 10 50 1350 + + 
4 8 a2 2040 + + 
5 6 18 2250 + + 
6 4 8 1728 = + 
7 2 2 686 + 
8 0 0 0 


Table 5-1 Example 5-3 for § = 0.5 
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Fig. 5-8a Stability and Dominant Root Regions of Example 5-3 
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Fig. 5-8b Root Configuration of Example 5-3 for 9 = 0.5 
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CHAPTER VI 
FEEDBACK COMPENSATOR DESIGN 
6-1 Introduction. 

The compensation by pure derivative feedback discussed in the last 
chapter usually adjusts only some of the coefficients of the character- 
istic equation. The effectiveness of this type of compensation is less 
than that if the compensator parameters enter all the coefficients. This 
has been shown by the reduced characteristic equation in Chapter ITI. 
Moreover, the tachometer feedback reduced KY and acceleration feedback re- 
duced K and so on; and a system can seldom be stabilized without adjust- 
ing the lower order coefficients of the characteristic equation. There- 
fore, if the compensator is chosen in the feedback path, the compensator 
often consists of a tachometer followed by a network. The tachometer 
effectively changes the coefficient of the first derivative term, while 
the poles and zeros of the network enter all the coefficients in the way 
of the coherent nature of algebraic equations. This chapter is a discus- 


sion about the design of such a compensator by the analytic method. 


6-2 Second Order Compensator. 
Consider Fig. 6-1 which is a second order system with feedback compen- 


sator consisting of a tachometer followed by a lead network. 





Fig. 6-1 Second Order System with Feedback Compensator. 
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The characteristic equation of the compensated system is 


S°+(84+ptKAa)S +(B+BPSt+tBp=O- 6-1 


where By = a, By = KK: In this characteristic equation, there are 


two variables in the coefficients, therefore two roots are arbitrary 
a K, and P have no restrictions. The reduced characteristic equation 
is of order one and has the form of S + ee = 0. The root coefficient re- 


lations are: 


Btpt+k,R, == 29Wnt+Co (6-2-1) 

B+ Bp == Wy t 26WnCo (6-2-2) 
2. 

Bp = Con (6-2-3) 


In equation (6-2), the dependent variables are p , k. and Co3 the in- 
dependent variables are f and Wn which are assumed to define the 
dominant roots. Co determines the stability and dominant root region, 
while p and Ra must assume values according to the dominant roots. 
Since equations (6-2) are linear functions of the dependent variables, 


it is convenient to write in the matrix form as follows: 


/ | = K> Re Se B, 
o B, -2¢0,| | P | |W -B, (6-3) 
| ee -Wy, & | 

: Oo | 


The three dependent variables can be solved readily, but only os must be 
expressed as an explicit function of the independent variables. By Cramer’s 


rule, one obtains 


C. 
B, _ Wy, 
° 28 Wy Bir 64) 
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e in this case is the actual real root. The stability and dominant 

root region are determined by equation (6-4) by varying < and W, , 
For design purposes, only those values of > and WJ, which are within 
the specifications are computed. When . and UJ), have been chosen, 


K,ka and P are calculated and the design is completed: 


Z 
Example 6-1: In the block diagram Fig. 6-1, K,K, = roe. aoa 123. 
Therefore B , = 10° and By = 23. From equation (6-4): 
(oe noe 
Co = 2x10° FW, 23 Wy, 
The stability criterion is Cc, >oO ine, 
> 2 
1O — W,, 
—————_—————- >O 
2X10*§ W,,-23W; 
or 


S Ps 5 2 
jO-W, >o and 2X10 SW, - 23W, — O 
These two inequalities define the stability region and is shown in Fig. 
6-2. If the dominant region is defined by Ci > 27Wp, then the 
criterion for this definition is 
lo wh, 
2X10 FW y,- 23wW* oa 
This region is also shown in Fig. 6-2. For g = 0.5, the range of W is 
e—4) =< 220. The values of C_, fP and K ka for the dominant roots in 
this region ( § = 0.5) are computed as shown in Table 6-1. The final 


decision of the design is to pick up a value which is favorable to all 


compensator parameters and other specifications. From Table 6-1 assume: 


WW), & p K, Re 
10 10° 10 9770 
20 5000 20 4960 
40 2480 38.4 2460 
60 1600 56 1580 
80 1200 77 1180 
100 922 92 900 
200 342 125 365 


Table 6-1 Values of P, Kka and C. 
of Example 5-1 for $ = 0.5 
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Stable root-region and Dominant root-region 
of Example 6-1 


110 





300 


200 


100 





The dominant roots are chosen to be § = 0.5, OG), = 100, then the third 
real root is -922, pP = 92 and K,ka = 900. The closed loop transfer func- 


tion is 


Oc oe 


7. ar 





oe tts) 
The closed loop pole-zero configuration of the above choice of the domin- 


ant roots is shown in Fig. 6-3, 


-922 


a 





Fig. 6-3 Closed Loop Pole-zero Configuration of 
Example 6-1 for the dominant roots 
i = 0.5, WwW, = 100 
From Fig. 6-3, it can be seen that the system is essentially second 
order because the real root is far away from the origin. The effect of 


te to the transient response is also easily visualized from the pole zero 


configuration. 


6-3 Third Order.Compensator. 

Fig. 6-4 is a block diagram of a third order system with feedback 
compensator. Assume the forward gain KK, is fixed from the steady state 
accuracy and the velocity constant must be greater than a given value Ko 
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Fig. 6-4 Third Order System with Feedback Compensator. 
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The characteristic equation of the compensated system is 
4 > 2 
S'+(B,tP)S +(8,+BptK,R,)5 +(Bptk, 2 48)5 +B p=o 6-5) 


where B= Cet. Pav 1Cilo , B= K,Kz 

In the above equation, there are three variables, namely R, . p and ay 
However, since K is specified, those three variables are not independent 
one to another, they are constrained by the specification Ke kK of the 
compensated system is 


K = Eee (6-6) 
i en ee 


Rearrange equation (6-6) 


__ PB. 
KR. e ak — PB, 


Substitute this expression into the coefficient of the first derivative 





term of equation (6-5), it becomes 





y a, B __ pBe 
BPtKR +B = Bp + e -PB+ 6 = - +t Bo 


Then equation (6-5) becomes 


Di 8. 
S*+(BtP)S + (Bt BP TK RIS t(—PtB)S+ Bp=o 7) 


The coefficients of characteristic equation now have only two variables, 
namely R, and Pp . The specification of K eliminates one of the vari- 
ables of the compensator. This is always true. The more specifications 
given, the less the freedom of the compensator. Since there are two vari- 
ables in the coefficients, two roots are arbitrary. The reduced character- 


istic equation is of order two and has the form: 
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Bt+p = 25W,+C, (6-8-1) 





ato. Pact: kK, Re aoa Ww Gee) On | (6-8-2) 

Bo e , 

KS P+ Bo == WC, + 28WnCo (6-8-3) 
2 

Bp — (a), ec. (6-8-4) 


where s and WW, are the arbitrary roots while a and Cy are the co- 
efficient of the reduced characteristic equation. As the four dependent 


) and Ry in the above equations are linear functions, 


variables C 


ir? Co 


it is convenient to express by matrix form as follows: 


zs WJ — 
2rPWn -l 8 | Co | 7B gy 
—W,, aestor Bk O p | | Bs 
O - We | O KR, | O 
ae = 
By Cramer's rule, one obtains 
eel one 
Bie na a (6-10) 


—— Bo 


o 2¢ B Wp ~ Kh _ ee 


The stability and dominant region are determined by C. and Ci: Choose 
dominant roots the parameters K Re F p and ZX of the compensator are 
calculated from equations (6-8) and (6-6) 

Example 6-2. For the system of Fig. 6-4 given a = 10, b = 30, KK, = 
21,000 K, = 60, damping ratio F of the dominant roots be approximately 
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s © a an 


n ¢ 
For the given values, Bo = 40 
By = 300 
B_ = 21,000 
Oo 


Substitute into equation (6-10) for f = 0.7 


2 3 
: 7 2.9 iB x [OW - 35x10 W, t coy 


For (), = 10, C= 1321 
From (6-8-4) P = 6.26 
From (6-8-1) C= B + P -26W, = 32.3 


Cc, and C, determine the secondary roots, they must be checked before 
going further. The reduced characteristic equation is 
5S + 3225 +1320 =O 
The roots of this equation are -16.1 + j 32.5. Compare the real part of 
these roots to that of the dominant roots ( Pw, = 7), it indicates the 
presumed dominant roots are dominant. 
From (6-8-2) 
KR, as i ae Ge Za a), C, - 8, TER =o 5 
and finally from (6-6) 


B= 226 , Z= 14,15 


The compensator becomes 
(323x%CSt (4,15) 
K, (St 6,26) 


The calculations above are for one point only. The choice of the 
dominant roots does not insure a satisfactory result. For more freedom 
and adjustment, the stability and dominant root region are calculated 


first as shown in Fig. 6-5. From Fig. 6-5, adjustments of the dominant 
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Fig. 6-5. Stability and Dominant Root-regions of Example 6-2 
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‘ Fig. 6-6. Closed loop pole-zero Configuration of Example 6-2 for 


§=0.7, W,= 10. 
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roots can be made, 


The closed loop transfer function is readily formulated as 





ae ms KK, (St P) 
Op Go > tO) (STS) 
The closed loop pole-zero configuration for 5 =e /, UW) = 10 


(dominant roots) is shown in Fig. (6-6). 

Example: 6-3. Given the same plant and the compensator as example 
6-2 except the specifications of KY is relaxed, Here are three free ad- 
justable parameters, three arbitrary roots are available. The control- 
led characteristic equation is the same as (6-5), the root-coefficient 


relation are as follows: 


eG, tS. t 5) (6-12-1) 
B+ BptkR, = Wy, +26 Wn(5,t54) + 3% (6-12-2) 
Bpt+kK,R,-tK = Wn (5,4 5,) + 2¢W, 5% (6-12-3) 
fae, S35, (-12-4) 


The four dependent variables are KK,» Pp >» Zz and Sys the three independent 


variables are § , W, and S,. Assume S, and S, are complex, and are 


nN 3 3 4 
4 
denoted by g and Ww eet UW)” is chosen as the dependent variable, then 


from (6-12-1) and (6-12-4) one obtains 
,¢ S K 
ie S 


Nn 


(6-13) 





oe) Wn — 





-B)=>0O 


This equation alone determines all the roots since it contains a Ca 


4 
9°’ , and Ww, . Assume § (dominant root) = 0.5, S ( secondary root) 
= 0.7, and substitute the numerical values of K = 21,000, By = 300 and 


= 40 into (6-13) and solve for Ge one obtains 


Wy’ — BtIXlo REGO 
n= Oe i+ 974x100 (WwW -40) (6-14) 
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Assume only 10 = W, = 20 are of interest, the results of compu- 
tations from Equations (6-12) and (6-14) are tabulated in Table 6-2. 


Notice in Table 6-3 when Gs 10 and UW = 23.52 is negative. 


10 Oh p KR, Z P/z 

10 270 348 62660 10.25 34 
2359 2.63 575 -7.25 -3.62 

15 109.2 128.5 9085 17.2 7.48 
21.4 

20 Ses 56.4 5324 223 2.53 
19.6 


/ 
Table 6-2 Example 6-3 for f= 0.5, § = 0.7 
For simple realization of the compensator, the roots of ‘J, = 15, 


wo = 109.2 are chosen. For this choice of the roots, the compensator 


is G = SCS an) 
e (S+ 128.5) 
where k = _Ke Re coe Ki = 1, k. = 0.433. tie four roots of 
t K, t 
the systemare (=0.5, Ww = 15, g = 007, t= 109.2. There 


is a zero at -128.5. From the closed loop pole-zero configuration, the 


roots i = 0.5, J, = 15 have high degree of dominance, 
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CHAPTER VII 
CASCADE COMPENSATION DESIGN 
7-1 Single Section: 
As cascade compensation is commonly used, the design procedure dis- 
cussed in Chapter IV, is applied. In Chapter I a third order system was 


analyzed, here a nth order system in the general case is considered. 





Fig. 7-1 Single Section Cascade Compensation for a 
nth Order System 


Consider the system of Fig. 7-1. Assume the characteristic equa- 


tion of the uncompensated system is 
n -| - “ = 
6 +B8,,5" +8,,5°°+ --- + BS +B8S+B=0 WD) 
For specified forward gain K and error constant (K Ko etc.) the 


three compensator parameters Z, p and C are not independent one 


to another, but related, namely 


aa 7-2 
oreo z 7-2) 


The characteristic equation of the compensated system is 


N+l : 
S tak B,tp)S" 4 cee ect ( B. + B.. p)S™+4 a + (Bx+ Bp)S (=>) 
+ Bee 2 


L19 





In equation (7-3), there are two variables PD and c& in the coef- 


ficients, then two roots are arbitrary. The reduced characteristic 


equation is of order (n-1l), and has the following form. 
n-l n-2 == 
S a Gate > , = >> +C,8 tC, =O 
If the two arbitrary roots are presumed to be the dominant roots of the 


compensated system and are expressed by f and OY, then the root- 


coefficient relations of equation (7-3) are: 
2 
BaP = Co, 
2 
Bo + BP = Ciw, + Co rf, 
z 
B, + Bp = O.4, + Cy2PW, + Co 


‘ é f ‘ 
t ; t 
’ ‘ : 4 


= alr B;P a Cw. t C204, + G (7-4) 


2 


B., + BP == Cane Ga iva + a 


na Nem! 


. a 


‘ i 
i ‘ 


; 
The N+l dependent variables ( p > xX ; Co Cyeee eC 5) are linear 


functions, then equation (7-4) can be expressed in matrix form 








a 
el O ae ee, ct ig | 2F Wy | C 
BS nt O O aa 2 Pw, on | Ww, | | ore 
Brnez O 2 fu, Ww Co O iP 
| i 2 | : 
B y-3 O j ’ “WW, : Ci | oO ; 
| 7 O | C2 — | : (7-5) 
! \ f ( ' | ea ' 
; ; | 
i i i } 
‘ai o ° | Ly | | 
, Oo ol 1 | | | e. 
B, O ey Pu, ’ ' \ 6, 
z ‘ { ‘ 
B, 5: ~2 Py “WwW \ t ; CO) 
‘ i 
B re ee en O Cae 
0 (0 n O O eu ees cee cl 





By Cramer's rule one obtains 


{ 


-BLRB + PB Wnt ZBWs +t --- FFB Wy 


(7-6) 
Bo o BW t fe b, Wis i we Bw, mis ve B., a, + 


where oe = —-'| 
ee 


ie [46° 
Y, = -4p48 Se 


i 


(727) 


{ 
a 
$= - [2599 - FJ 

J 

The other dependent variables can be expressed in the same way but this 

is not necessary as far as computational purposes are concerned. When Co 
has been calculated, the other variables are readily evaluated from equa- 
tions (7-4). Choose y , evaluate the stability and dominant root 
interval of Wry . Within the dominant root region, choose (J, and evalu- 
ate p and Q& . If the network is realizable, the design wrok is 


completed. 


oO 
Example 7-1. Given G = oer and K is not to be reduced. Then 


the characteristic equation of the uncompensated system is 


YY 
+ 
> 
Oo 
© 
HI 


0 

B. = 1, B, = QO, B , = 400, B, = O where i = 2 

The order of the compensated system is three and that of the reduced 
characteristic equation is one. Then 


C, =1l, C. =0Q0 wheni => 1 
1 i 
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From Equation 7-4: 


- 


ao CoWy 
BX + 8,p = nee + C,29W, 
ar B, Pp — C,25Wn + Co 


Substitute the values of B's into the above equations, obtains 


4G.00p = Cw 
400k = Wy +2PWnCo 
ee CO ee 
From equation (7-6) and (7-7) 
6, ( 8, - 2f Wn) 
Bea 


C, = 


Substitute the values of B's, obtain 


Wn 
ee 


WJ, — 400 
For stability the criterion is (C,»> 0, namely 
GOO FW 


2 
u,? — 4.00 


The numerator never becomes negative, then the inequality becomes 


2 
W, => 400 


(7-8) 


(7-9) 


(7-10) 


(7-11) 


GD: 


Define the dominant roots such that the real part of the secondary roots 


must be at least twice that of the dominant roots, then 


§00 Fn 
@,, > 25W, , “w= 400 


Then the criterion for dominant roots is 


doo< W,, < S800 


—=> 2fW), 


Those regions are shown in Fig. 7-2. 
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Within the dominant root region, for . = 0.5, and r = 0.6, the 
compensator parameters are computed from equation (7-8) for different 
values of Wy and tabulated in Table 7-1. Table 7-1 includes not only 
the compensator parameters, but also all the roots of the system for dif- 
ferent choice of the compensator. Assume the dominant roots are chosen 


as y = 0.6, WW, = 30, then from Table 7-1, C. = 28.8, P = 64.8, 


441( 84148) 
KX = 4.41, = = = 14.8, The compensator is G_= ae 
aos 
Co ie n Co ie A 
220 242 ee 264 2911 14 
54.6 78.5 4.72 67.7 94,3 5.34 
37.7 63.4 4.12 45.3 76.1 4.6 
29.2 Dine 4 Bak 68.5 4.41 
24 54 4.05 26me 64.8 4.41 
20.5 54.6 4.2 24.6 65.6 4.53 
L933 Spl) 4.53 oS ye 66.8 4.87 





Table 7-1 Computation Result of Example 7-1 
The closed loop-pole-zero configuration for this choice of the dominant 
roots is shown in Fig. 7-3(a). In Fig. 7-3(a) it can be visualized that 
the peaking time is small but the overshoot may be too large because the 
zero is too close to the arieine., If < = 0.6, “, = 24, are chosen as 


the dominant roots, then from Table 7-1, oe 67.7, P = 94.3 ,& = 5,34 


4 


and 4 = 17.7. The compensator is 


634( $+iht) 
oa ($4943) 


The pole-zero configuration for this choice of the dominant roots is shown 
in Fig. 7-3(b). For this choice of the dominant roots the peaking time is 
still small, but the overshoot is reduced. 
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Fig. 1-2 
Stable and Dominant Region of Example 7-I 





Fig. 7-3(a) 
Closed Loop Pole-Zzero Configuration for Dominar* 
Roots Si somo. WW), = 30 for Example 7-1. 


~~, 





— 


Rit 


Bee f—5\D) 


Closed Loop Pole-zero Configuration for Dominant 
Roots g = 0 on W), = 24 of Example 7-1. 


tn 
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we ea 
7 5$€$+2)(5+5) 


reduced. The characteristic equation of the uncompensated system is 


Example 7-2: Given G and KY is not to be 


so 4 1B + 102 + 300 5 0 
From Equation 7-4: 
iia = Cre 26 Wy 
IO+ TP = W, tle +2G WHC, (7-12) 


Z 
JOP + 300A = CW, +2°WnCo | 300P = Cow 
From Equation (7-6) and (7-7) 


> -300 [10 +14 FW (1-49 Wr | (7-13) 
a B20 - (W425 Zoe 
For i = 10). 3, 
— ( T&, ~10) BOO 
° Tue = Wd, - 300 
For =. 9 


C= 300 (0.AW), + 4.2%n ~10) 
oT TW, —2.6w), - 300 


To determine the stability region, the criterion is 


7-14 
eno and eee a 
Define the dominant roots by a = fW),,, from Table 2-7 of Chapter IT, 


the transformed coefficients of the reduced characteristic equation are: 


1 C,-2Q@ 


i 


Deen eres 


Oo 


The criterion for dominant roots by this definition of dominant roots is: 


Digese and Pee 5e (7-15) 
Computations for a 0.3 and 5 = 0.5 from Equations (7-13) and (7-12) 


are shown in Table 7-2 and Table 7-3 respectively. 
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Wn ee p C, x D, dD, 
0.2 9.13 0.001217 7 0.0046 + + 
0.4 8.22 0.0044 6.88 0.0101 + 1 
0.6 7.25 0.0084 6.81 0.0166 + + 
0.8 6.34 0.0135 6.49 0.0235 +r + 
r0 Bre o 0.0176 6.41 0.0313 + + 
2 4.16 0.02 Saas 0.0397 + r 
1.4 2.97 0.0193 6.21 0.049 = r 
TO re7l 0.0146 6.07 0.057 + + 
ies 0.397 0.0043 8.86 0.097 + + 
1.9 -0.312 -0.0038 a 13 + - 


Table 7-2 o = 0.3 for Example 7-2 


Wr, ee : Pp ej oh D, jon 
0.6 5.84 0.007 6.4 0.019 + +r 
0.8 4.46 0.0095 ba? 0.0191 + + 
1.0 3.07 0.01022 6.01 0.0216 + + 
ea 224 0.00904 6.02 0.0300 + - 
1.2 ie 30 0.00748 6.05 0.0322 + - 
Le 0.87 0.0049 5.75 0.036 r - 
1.4 0.1925 0.00125 5.61 0.0375 + - 
ae, -0.48 -0.0036 5.46 0.0383 + - 


Table 7-3 P= 0.5 for Example 7-2 
From the criteria (7-14) and (7-15), the stability and dominant region 


for F = 0.3 and § = 0.5 are determined from the computation tables. 


Stability Dominant ( a= fw, ) 
F = 0.3 U),< 1.8 Ge las 


Since this computation is a progressing process, when the stability 
and dominant root region have been located, the compensator parameters 
also have been calculated already. The computation tables not only in- 
clude the compensator parameters, but also the secondary roots. The 
coefficients (D's) of the transformed equation are shown signs only in 
the compensation tables, because in this case only signs are needed to 


determing the dominant region. 
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Assume the dominant roots are chosen as § = 0.3, W), = 1.0, then 


from Table 7-2, Co = 5.28, C, = 6.41, P = 0.0176, @ = 0.0313 and 


1 
—— z = 0.563. The reduced characteristic equation is 


$7 + 6.418 + 5.28 = 0 


Then the secondary roots are S = -5.44, -0.97. The closed loop pole- 
zero configuration is shown in Fig. 7-4(a). In Fig. 7-4(a), all the 
secondary roots are to the left of the dominant roots. 

Consider Table 7-2 and 7-3, they indicate OQ — 1, this implies the 
compensator is a lag network. For a lag network the zero is close to 
the origin. Then a root to the right of the dominant roots is allowed. 
This dipole near the origin has negligible effect on the transient response 
if the root is close enough to the zero. In this sense, the definition of 
the dominant roots may be adjusted as such that all the secondary roots 
must be to the left of the dominant roots except the root which is close 
to the zero. 

In this example: if 9 = 0.5, W, = 1.2 is chosen as the dominant 
roots, from Table 7-3, C, = 1.56, C, = 7.25 pb = 0.00748, & = 0.035 


Pp 


and x= oe 0.214. The reduced characteristic equation is 


5 + 7.258 + 1.56 = 0 
Then the secondary roots are -7.03 and -0.225. The closed loop pole- 
zero configuration is shown in Fig. 7-4(b). Compare the pole-zero con- 
figurations of Fig. 7-4(a) and Fig. 7-4(b), it can be easily seen that 
the scheme of Fig. 7-4(b) has better performance than that of scheme 


Fig. 7-4(a). 
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Fig. 7-4(a) Closed Loop Pole-zero Configuration for the Dominant 


aN 


ss 
~ 7.03 


Roots P = 0.3, WwW, = 1 of 


Example 7-2 
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Fig. 7-4(b) Closed Loop Pole-zero Configuration for a Dipole | 


Near the Origin of Example 7-2. Dominant Roots 


Wy = a2: 


1°70 





7-2 Stability Boundary Lines. 

It has been pointed out that the stability boundary lines for a 
cascade compensator are the same boundary lines of relocation zones in 
Reference 7. In the single section zone if the compensator is to be 
realized by an R-C circuit with pole-zero on negative real axis, then 
the necessary and sufficient condition in all the roots are within that 


: . x(S+ 2) 
zone. That is in the single section zone the compensator, Bape ee 


c Cr 
always can be realized by an R-C circuit for any choice of the wots oe 
viding those roots are within that zone. In the multiple zone, the roots 
being in that zone is a necessary condition but not a sufficient condi- 
tion. 

For higher order systems there may be several zones for a single 
section compensator. But for dominance, only the zone which is the 
closest to the jyW axis is of interest. It has been proved that putting 
a pair of complex roots within the lag area and close to the stability 
boundary lines (uncompensated root-loci) implies a closed loop dipole 
close to the origin. By the same reason, putting a pair of complex roots 
within the lead area and close to the boundary line (pseudo-root loci) 
implies all secondary roots are far away from the dominant roots. The 
question of how far away and how close they are, the dominant root bound- 


ary lines of the definition a = GF Wn serves a guide for both condi- 


tions. 


7-3 Double Section Cascade Compensator. 

If the dominant roots are expected to be within the multiple sec- 
tion zone, at least two sections are necessary. In general, the design 
is carried out by two steps. First choose an intermediate dominant root 


which is close to the stability boundary lines. Form the transfer function 
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obtained by the intermediate dominant roots, then repeat the process. 
The working formulas for the open loop transfer function with one zero 
in it have been derived in Appendix 2 for an original nth order system. 
The design by intermediate dominant roots has two advantages. First the 
quadratic term as discussed in Chapter 3, equations (3-3) is avoided. 
Secondly the compensator can be designed as a double lead section, a 
double lag section or a lead-lag network. The intermediate dominant 
roots should be chosen close to the final dominant roots because the 
regions which are close to the dividing lines (both uncompensated root- 
loci and the pseudo root-loci) are dominant root regions. 

4.00 


Example 7-3: Given CS ee 


on a pair of dominant roots | =0.7, Wy = 20, and K, is not to be re- 


, the system is expected to operéte 


duced. 

From Example 7-1 and Fig. 7-2, WW, = 20 is on the stability boundary 
line, then double section is needed. From Table 7-1, choose the inter- 
mediate dominant roots P= O24 W = 22, which is close to the final 


dominant roots ’ = 0.7, W = 20. The transfer function becomes 


52800 St18S) (7-16) 
s*C S$ + 242) 


G= 


From Appendix 2 equation A-2-4, the four equations are 


P-—C, — 259W,- B> (7-17-1) 
2 So ait 

Bp + 220 - Co -2 fn, = W,-8, +, (7-17-2) 
i =- = 

Bp + Bak - 26 yo —-W,C, =O (7-17-3) 

(7=17-4) 


2 
BP - Wi Co =O 
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From A-2-6: 


B. | B,- 25 6,9, = (1-4p Jun F = [8,- Bw), trp, | 


a | B.- Bw, top, | ir = eee Bo “Ly, | 


(7-18) 


c 


where Bo? B, and Bo are the coefficients of the characteristic equation 


1 
of the open loop transfer function (7-16), (-C) is the open loop zero. 
The stability boundary lines are shown in Fig. 7-5a by the rules of 

reference 7. From Fig. 7-5, it indicates that the roots of , = 0.7, 


U),, = 20, are within the single lead section zone and close to the pseudo 


root-loci., Then they are expected to be dominant. Substitute C = 185, B A = 


97500, B, = 5280, B, = 242, : = 0.7, W,= 20 into equation (7-18), one 
obtains 
C_ = 6340 
Oo 
From (7-17-4) p = 26 
From (7-17-1) Cc, = 240 
From (7-17-3) AO = 1.4 
Then Z = leo 


The secondary roots are determined from the equation 


6*+ 2405 + 6340 =O Cae 


The solution of (7-19) is -210, -30. 


Z $4/3.6 
The compensator for the second stage is GP == misses Cae va 
(S+2G) 


The overall compensator is 
— 18.5 CS + 18.6) 

G. = (S + 242) CS +26) 
The closed loop pole-zero configuration is shown in Fig. 7-5b. From 


Fig. 7-5b it can be visualized that the response is speeded up by the 


two zeros while the overshoot is small because of the approximately equal 
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Fig. 7-5(b) 


Fig. 7-5(a) 
Relocation Zones of 


G = 62§O ( § +185) 
6? ( $4242) 
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Final Closed Loop Pole-zero Configuration of 
Example 7-3 
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distances of the zeros and the dominant roots to the origin. 


7-4 Mixed Compensators. 
For a single section cascade or feedback compensator it has been 
shown that the adjustable parameters are linearly dependent (Matrices 6-9 


and 7-5), In general a compensator has the transfer function 


nie ears a: =) aie 
i “oe. | —— 
C OS + oan ees 
The coefficients, AQ,» A, » -- Gp and Deo bg are linearly 


dependent for either a cascade or a feedback path. The solutions given 

a set of roots are always real numbers. If there are more than one compen- 
sator and in different loops, the adjustable parameters are not linearly 
dependent. The solutions for a given set of roots may be complex numbers 
and consequently the compensators are not physically realizable. However, 
the subset of roots which gives the solution of real numbers can be found 
if such a subset exists in the stable root region. Note the coefficients 
of the reduced characteristic equations are always real numbers for any 
choice of the roots and for any type of compensators. The reason for this 
fact is due to the complex conjugate roots, Because of this fact the stable 
root region and dominant root region are always defined if they exist. 

The physically adjustable parameters also appear as real numbers in some 
transformations for any choice of the roots, but the individual parameters 
may not be real numbers. Consider Fig. 7-6 which consists of a single 
section compensator in cascade and another in the feedback path, Cr is 

the given plant and Z, , Z,, PP, » f and Ri are adjustable parameters. 


From a signal flow graph, the characteristic equation is 


(StP)(StP)Q +R StE)(StB)P + RS (St 2.)(5+ P) P=O 
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(7-21) 












K.¢$4+2,) 
(5+ P) 





£5 (St 22) 
(S+PB) 


Fig. 7-6 Mixed Compensators 


O(St2) 1” 
Y (St P) Y SCS +5) (S+10)* 


RS 


Fig.7-6(a) Mixed Compensation of Example 7-4 


134 





Consider the first term of (7-21), it can be expressed as follows: 


(S4+P)(StRJIQ = [St CPtB)s + PRI @ (7-22) 

Let 
et =O (7-23) 
PRP = oa : (7-24) 


where 6 and Q are always real numbers for any choice of roots of (7-21). 
However, , and Pp, are not necessarily real. To find the real solution 
of p, and BR > the roots must be subject to some restrictions. In this 
case it is 
2 
b >4Q (7-25) 
The inequality (7-25) defines the realizable region. For a passive net- 


work, additional constraints are necessary, that is 


iO Gio (7-26) 


The same constraints of the roots can be applied to other cases. 





Ss 
Example 7-4. Given G — 10 , The velocity error 
S$CStS)CStHIO) 
constant Ky is not to be reduced. The band width is less than 20 rad/sec. 
From the specifications dominant roots of ¢ = 0.5 Wz 10 are chosen. 
For cascade compensation alone, more than two lead sections are required 
by the Ross Warren method: A scheme of Fig. 7-6(a) is attempted. Be- 


cause of the gain requirement, only second derivative feedback is used 


and the parameters & , zZ and fp are related by the following relation 


wee (7-27) 
a 


There are three free parameters (a, pos an ), three roots are artit- 


rary. The constant controlled characteristic equation is 


oo 4 (264+) 5+ + (2004 500 P +10°Ka)S + ($00 +t 200P tO Kp)s” 
4 (scoop tio’) 5 + 10 p =O (7-28) 


> 





Let (-S,); (-S,) and (-S,) be the arbitrary roots and define the 


following terms for simplification. 


eet. tS, L7- 29-1) 
peo ee, +S; 5, (7-29-2) 
F; = S513 (7-29-3) 


The root-coefficient relations are: 


S00 +P = F,4C, (7-30-1) 
HOO +iGOOP 410k. =MCs + CF, + F, (7-30-2) 
Be coop top, Pp = Coke Fs (7-30- 3) 
25Pp +hio’ = Col, + CFs (7-30-4) 
jo P = C.F, ve 


Where Cy and a are the coefficients of the reduced characteristic equation 


which defines the constrained roots (-S,) and (-S_) by the equation 
2 
et CS Ge =e Gaal) 
From equations (7-30-11), (7-30-2), (7-30-3) and (7-30-5), er is 


obtained: 
. _ 2 io oes oe lg 
F, (FF -25F +10) C, — 10 [FE -2665 +10 Ele 


tio’ [ EF, - 256 -F, +500] =O a 
Assume (-S)) and (-S,) are the dominant roots and denoted by ° and GW, , 


(-S,) is arbitrary and must be greater than (W, . Assume S, = 2f Wn 


For ° =o LW) = 10 from (7-29) 


F, = 20 
F, = 200 
F, = 1000 
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ee 9 errr @ 


— a. a — mete & 


| cleanly 





Substitute into (7-32) one obtains, 


Cee Z050R 31.6 
Take Co = 2080 
From (7-30-5) P = 20.8 
From (7-30-1) C, = 2900 
From (7-30-2) 10 Ry = 2076 
From (7-30-4) OX = 4,32 
The reduced characteristic equation is 


The 


The 


The 


The 


5° + 2685 + 2080 =O 


constrained roots are “Sy “5. —eale.9 +eme) 43.5 
five roots are: 

“Sy> “5. = -5 + | 8.66 

“5, = -10 

“S,> “Se = -12.9 + j 43.5 


compensator is 


— 432 ($4483) 
Te $+208 


a O02 Ole 


closed loop pole-zero configuration is shown in Fig. 7-7. For other 


choices of the arbitrary roots, better results can be obtained. 


X 


30 





Fig. 7-7 Closed Loop Pole-zero Configuration 
of Example 7-4, 


Ly 





VIII. CONCLUSIONS 
8-1 Summary of Results. 

A new method for design of control systems in the time domain has 
been described. The method is applicable to any system configuration and 
is able to handle any number of variable parameters. The design procedure 
is based upon a set of simultaneous algebraic equations, and is well suited 
to implementation by digital computer. The method is flexible and has more 
freedom in design than other methods. The values of adjustable parameters 
are computed from analytical expressions, so any desired degree of accuracy 
can be obtained. 

The main advantage of this method over those methods available in the 
literature lies in the fact that the method is a direct syntheses from the 
closed loop pole-zero configuration. No trial and error is needed. The 
specifications and restrictions of the compensator confine the arbitrary 
roots within a region on the S-plane, once this region or regions have been 
determined, the rest of the design work is just a computation. The compen- 
sator obtained by this computation guarantee that the specifications and re- 
strictions are satisfied. 

The main difference between this method and other available methods 
is that the domain of the system can be any number of variables. Theoreti- 
cally, the more variables available, the more design freedom one has. In 
the root-locus method, only one parameter is varied others must be fixed. 
In Mitrovic's method, only two parameters are varied. This method is able 


to handle any number of variable parameters simultaneously. 
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8-2 General Description of the Method. 

The design procedure can be summarized as shown in Fig. 8-l. The 
analytical expressions of the dependent variables of commonly used com- 
pensators have been derived in the previous chapters. [In general for the 
case in which coefficients of the characteristic equations as variables, 
the dependent variables are always in linear relations one to another. 

A compensator of the transfer function as (8-1) is in the forward path or 


in the feedback path alone, tne coefficients a's and b's are always linear 


G= pei ek aa tO An ST bee te (8-1) 
e 5. ae SiR cee. De 

relations in the partitioned formulas. Those linear algebraic equations 
are always consistant and therefore always have a solution. For the case 
in which compensators are in both forward and feedback paths, a quadratic 
expression is always obtained. But this quadratic expression does not 
imply the dominant and realizable regions are smaller. It shows only the 
natural algebraic relations. For example, the coefficients of the reduc- 


ed characteristic are always real no matter how complicated their functions 


are. 


8-3 Extensions. 

The method can be extended to sampled data systems and adaptive sys- 
tems. The modifications are the different interpretations of the variables, 
For example, in the case of sampled data systems, if A transform is used 
then the criteria of stability, dominance and realizability are different. 
However, the general idea of approach can be applied. 

In the previous chapters, the arbitrary variables are chosen as the 


roots of the system. This choice is not necessary. Theoretically, any 
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combination of "r'' variables can be chosen as the arbitrary variables such 
as a combination of roots and variable parameters. 

This method can be applied to any mode of linear operation in dis- 
continuous operation. By proper choice of the compensator, the dominant 


roots of two modes can be chosen in a relation as the designer wants. 
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APPENDIX I 
CONSISTENCY & LINEAR DEPENDENCE OF DERIVATIVE FEEDBACK 
I. Introduction. 

The coefficients of the characteristic equation of a linear control 
system are functions of the system parameters. If all the coefficients 
can be adjusted by a set of the parameters, then the roots of the system 
can be controlled. Conversely, if the roots are chosen, then the system 
parameters are forced to have certain values. The conditions of using 
system parameters to satisfy the coefficients of a cnavacte tenet equation 
for a given set of roots depend upon the number of adjustable parameters 
and the consistency of the simultaneous equations. Theoretically, if all 
the derivatives of a system can be obtained from the output of the system, 
then by feedback as shown in Fig. A-l-l, the system roots can be chosen 
arbitrarily providing no restrictions are placed on the forward gain K 
and the control parameters h's. This can be shown as follows. The charact- 


eristic equation of the system with feedback is: 


Be (Gh +a,.)57 "4 Lt Kh tals + KK =O (1-1) 


Assume the roots of the system are chosen from the considerations of 
the specifications, and the required characteristic equation is formulated 


as equation (A-1-2) 


Seen Ss ee SS ars Pas HA =O e(a-1-2) 


For the system to have the chosen roots, then the corresponding coeffi- 


cients of equations (A-1-1) and (A-1-2) must be equal. Then 
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h,..> + Ans + aie a + 4,5 


~~ 


\ 


ae 
Fig. A-l-1 Block Diagram by Feedback of Derivatives . Cr is 
Transfer Function of the Plant . K, and K, 
are the System Gains. 


143 





KK, = H » 

ore (A-1-3) 

K, hn na t Qn — A n-' 

Equation (A-1-3) always has a unique solution if there are no re- 

strictions on K's and h's. And in the particular case, each parameter 
h controls one coefficient independently and the forward gain controls 
the constant co-efficient. If the control parameters have certain restric- 
tions, then the system roots cannot be arbitrarily chosen, but have to be 
subject to these restrictions. For example, if Ky Ko» is specified from 
the consideration of steady state accuracy, then the product of the roots 
= 0), then the sum 


must equal to K If h cannot be obtained (h 


82° = 


of roots must be equal to aa 


1 


1° The restrictions on the control para- 
meters restrict the roots in a certain constraint. The more restrictions, 
the higher is the degree of the constraints of the roots. Ina practical 
system, the higher derivatives are very noisy and very difficult to obtain. 
This is a common restriction for derivative feedback. 

In Fig. A-l-l, it was assumed that there were no zeros in the forward 
path. If one zero is in the forward path, then the allowed highest order 
of the feedback derivative is (n-2). In general, if ‘'yn'' zeros are in the 
forward path, only the derivatives of the order up to n-(M+l) are allowed 
to feedback. This can be proved as follows 


Refer to Fig. A-1-1, assume m zeros in C7 _, and the highest deriva- 


tive has an order of ,(  . Then the characteristic equation is 


Seca, thy, |S 2a eo (A-1-4) 


144 





The subscript ( tm ) of ao indicates the highest order of the feed- 


back derivative and the number of zeros. From equation (A-1-4) 
r+im =z n—-] 


r< N-Cmtt) 

Therefore if m zeros are in the forward path, then only n-(in+i ) 
control parameters can be introduced, if the feedback path is as shown in 
Fig. A-l-1l. However, other schemes of lower order derivative feedback are 
possible to introduce (n-1) control parameters and consequently control all 
the coefficients of the characteristic equation. The following sections 
are discussions about these possibilities. In the following sections, two 
terminologies are defined as follows: 

Control variable is the variable of the system which can be measured. 


Control parameter is the adjustable parameter in the system, 
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Zs Feedback of the First and Second Derivative fromthe Output to 
Several Points in the Forward Path. 

Since the first and the second derivative are usually allowed to 
feedback for high order systems and can be obtained, only these two deriva- 
tives are considered, The object of this discussion is to adjust all the 
coefficients of a nth order characteristic equation to desired values. In 
order to accomplish this, n control parameters are needed. As the forward 
gain K is considered as one control parameter, (n-l) control parameters 
from the feedback path are necessary for this particular scheme, Using 
first and second derivative feedback around the plant intorduces only two 
control parameters. However, if there exist junction points between the 
energy-Storage elements in the forward path, such that signals can be fed 
in or picked up, then additional control parameters can be introduced by 
feeding derivatives to those points. As two control parameters have been 
introduced by feeding derivatives to the error channel, then the neces- 


sary number of points between the energy storage elements in the for- 


n-3 
2 


derivatives can be fed into each node. Briefly, this can be summarized in 


n-2 





ward path is (if n is cdd and if n is even) since two 


a 


the following two statements. 
(1)* If a nth order system has at least nes available points 
2 


(if n is odd, and iS if n is even) between the energy storage elements 
in the forward path, then all the coefficients of the characteristic equa- 
tion can be adjusted to arbitrary values by introducing first and second 
derivatives from the output to each point providing the two derivatives 


are defined and the feedback paths are proper. 


*all statements in this discussion may be considered as rules. All of them 
are necessary conditions, but not sufficient conditions. 
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(2) If a nth order system has at least n-2 available points 
between the energy-storage elements in the forward path, then all the co- 
efficients of the characteristic equation can be adjusted to arbitrary 
values by introducing only first derivative feedback from the output to 


each point, providing the feedback paths are proper. 
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Di Formulation of Characteristic Equation for Pure Derivative Feed- 


back. 


From the signal flow graph, the determinant of a system can be 


obtained by the following rule. 


(«) = (2) 
Ae je” +2" 
nr) eS K 
where i = product of the kth possible combination of “2 non- 
K 


touching loops 
[, = loop transmission 
For this particular derivative feedback system, two terminologies 
are defined as follows: 
Loop gain = numerator of each loop transmission. 
In Fig. A-1-2. Loop gain of Ly = h,P, 
Loop gain of Ly = hy oP 
Loop gain of Ly = haP PoP, 
Loop gain of Ly = PLP oP, 
Characteristic function = denominator of each of the forward transfer 


functions, or the product of them. 


Consider Fig. A-1-2: from the graph rule 
ire it i, + 2h bs, tee 


Ph PF, he Pte ae Pelee 
sa + 3 23 p22 
TQ" aa ea * aaa 


| 


\ 


QQ,Q; T Q,Q, (Ph) + @mcR eh, } 
slemteetcenoe te ora kes (A-1-5) 


The first term G), QQ, + last term P,P,P, = characteristic equation 


of the original system, 
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Fig. A-1-2 Formulation of the Determinant 


P.(s), Q;(s)_ h.(s) are Polynominal in S$ 
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The second term = loop gain L, X nontouching chara. function. 


1 


The third term = loop gain L, X nontouching chara. function. 


2 


The fourth term = loop gain L. X nontouching chara. function. 


3 
From equation (A-1-5), the formulation of the characteristic equa- 


tion for pure derivative feedback may be stated in the following rule. 


The characteristic equation of a pure derivative feedback system = 


(\) 


characteristic equation of the original system + = e loop gain X 
ei : : (2) ; 
characteristic function of nontouching loop + =. e loop gain X 


characteristic function on nontouching loop. 

Example A-1l-1: Devise a derivative feedback scheme of Fig. A-1-3, 
such that all the coefficients of the characteristic equation can be ad- 
justed. Assume the required characteristic equation from the specifica- 
tions is 

So 4 As tA,S + Ars*+AS+tA,=O 
find the relations of the control parameters and the coefficients of the 
given characteristic equation. 

Since intermediate points are available, then n-2 = 5-2 = 3, From 
rule (2), it is not possible to use first derivative only. 

By rule (1) a = =? = 1. Then a combination of first and second 
derivative can meet the requirement. There are many possible ways to 
accomplish this since only 4 control parameters are needed. 


Two schemes are listed as follows: 


2 
1 1 9 h,s + h,s Hy a h,s 


it 
ae 
wn 
jee) 
it 


(a) H 


2 2 
1 1 , hos H, = Bs + hs 


il 
a 
n 
ei 
it 


(b) H 


Assume scheme (a) is chosen. The chara. equation is formulated in 


a tabulated form for convenience. 
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Fig. A-1-3. Example A-1-1 - Mark "X" denotes signals are allowed 
to feed in or pick up 





St 
Required Chara. Equation S| 
Original Chara. Equation 

Gain |_, x Chara. Nontouching 


2 f 
h, 5 Ki K, Ky (ste) 


Gain L, 


(h,5 + h,S)K, Ks *CS+ Pi) 


x Chara. Nontouching 


4 3 2 | 
(Aas | A3S As 


i 


| POA ae 


Za Palate ls oe | ie A KKK, K3 
+ KKK 
| 


KK, K; Na KK, KX h4 | 


Kaha | KP he Ik keh, 


Gain Es x Chara. Func. Nontor, | K hy K3h,(p +2, Kk (pe + gree, 


hs Ks x (STP GER) SP) 


sate) BPs+hP )h, 


The equations for the coefficients are 


K, ri 2 = Ag 


Za PC; Ts K, Ka ha ie KL h, te k, (pt Pt Ps)h,= Az 


= aighe t KK, Ksah, + K, Ks BAe +k, ks hs t lene eis tig ie 
Deir iicick, t KK PA, +K,PRRGh, = Ay 


Bae iiciks =A, 


The first four equations of (A-1-6) can be written in matrix form: 


K; o 
Ks CF +t, +R) KK; : 
K3 ( i (2+ pie PP.) K, K; P KK 
K 3 fies O KK, 


IS 


2 h, Ai = Eee 
KK, K3 | | he A,-2f,f- 
KKK] | Neg A,- =P 





In (A-1-7) whenever the determinant is not zero, there always exists 
a solution of the control parameters h's for arbitrary values of A. 

In the above example, the table of the characteristic equations gives 
a clear picture of how the coefficients are controlled by h-parameters. 

Also equation (A-1-7) shows that the h-parameters are linearly dependent 
to the coefficients of the given characteristic equation. It is this 
linear dependence which makes the arbitrary adjustment possible. In 
equation (A-1-6), the constant term is controlled by the forward gain only. 
This is the nature of derivative feedback, because the derivatives control 
the coefficient of the terms of a differential equation. Any derivative 
feedback from the output has the scheme of properties of linear dependence 
of h-parameters and the independence of forward gain K. 

Example (A-1-2). Fig. A-1-4 is a third order system, one point between 
the energy-storage elements in the forward path is allowed to feed in signal, 
then the system roots may be arbitrarily chosen. Assume the roots to be 
-3 and -0.5 + j l. K> h> h, can be determined. The characteristic equa- 
tion for the given roots is; 

S + 4S” + 4.258 + 3.75 = 0 

Main loop 1 1 0 K 


i 


L,: loop gain = K) hos 
chara. function 


not enclosed = 1 Kh, 


L,: loop gain = h,s 
Chara. function 
not enclosed = s + l h h 
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Fig. A-I1-4 Example A-I-2 
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Fig. A-I-5 Example A-I~3 
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Set the corresponding coefficients equal, variable parameters are solved. 


Ky = 3675, h, = 3, hy = 0.334 


Example A-1-3. In Fig. A-1-5, it is assumed that signal can be fed in at 
the point "X" in the forward path, and the forward gain of the system is 
1500 based on the steady state specification. According to rule l, ws = 
then by first and second derivative feedback, the system roots can be chosen 
arbitrarily providing the product of the roots is equal to 1500. 

The uncompensated system is unstable. The feedback paths are chosen as 
shown, and the roots are assumed to be S = (-3 + j5), -6, and -7.36. 

The characteristic equation for the chosen roots is 
s* 4 19.365° + 158.5° + 7208 + 1500 = 0 


The system coefficients 


Uncompensated system = l 10 St 30 1500 
i h)sth,s ) 1500 1500h, 1500h, 
L,: h; S 30 ($+3)(S+$) 30h, 240h., 450h., 


The coefficients relations are: 


30h, + 10 = 19.36 


1500h,, + 240h,, + 31 = 158 

450h., + 1500h, + 30 = 720 

Solve equation (A-1-8), get (A-1-8) 
hy = 0.367 
h, = 0.0347 
h, = 0.312 
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Di Feedback from Intermediate Points. 

It can be shown that if the output of each block in the diagram (A-1-6) 
is chosen as the state variable of the system, then there exists a linear 
transforamtion between the state space and the phaSe space (output and all 
its derivatives) of a system. In Fig. A-1-6, the state variables are 


chosen as shown, the relationship between state space and phase space is 


Ka Xx Kb Xo f Kiaioe 
CC SAG Stb -, 


Fig. A-1-6 State Variables 


derived as follows: 





X, =o (21-91) 
a | _ ve cl A-1-9-2 
se = ey X2= xX, = t ( ) 
cca _ = + (48 (A-1-9-3) 
emote ete Gs = Kola + «(Ge 2) 


Equations (A-1-9) are the required relations. It shows that the out- 
put of each block is a linear combination of the output © and its de- 


rivatives. 
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With this concept, therefore, the feedback from intermediate points 
is equivalent to derivative feedback from the output. Fig. A-l-7(a) is 
equivalent to Fig. A-1-7(b). 

This equivalence also can be shown by block diagram manipulations. 
Because of this equivalence, if the signals flowing in the intermediate 
points can be measured, then those signals can be taken as control vari- 
ables and the analysis is the same as for derivative feedback from the 
output. Moreover, if the derivatives of those signals are defined and 
can be obtained, they also can be taken as control variables. 

However, if consistency and linear dependence are considered, the 
feedback paths and the number of control parameters have many restrictions. 

(1) Consider the scheme of Fig. A-1-8, it introduces sufficient 
control parameters but the coefficients of the characteristic equation are 
not linear combinations of the h-parameters. It consists of two indepen- 
dent nontouching loops, the contribution of these two loops to the whole 
system is the product of the two loop transmissions. So the coefficients 
are quadratic functions of the h- parameters. This can be shown by formu- 


lating the characteristic equation as follows: 


5° 57 5! oe 
Original system 1 OAt+b ab Ka kK, 
L, : k,h,$ (5+a) Kh, Kh, 
L, Kh, (sb) Kh, Kb hz 
L, x L, KK hh, 


If the required characteristic equation is 


5° ae A 5? + ALS +A = (0, 
Z 1 oO 
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Fig. A-1-7(a) Feedback from Intermediate points. 
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Fig. A-1-7(b) Equivalenence of Fig. A-1-7(a) 
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Fig. A-1-8 Quardrtic Function of h parameters 






| Kb 
S(S+ b) 


Fig. A-1-9 Example A-1-4 
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then 


Kh, + Kah, = A,- (at b) 


Kaoh + Kb h, + KK hh = A-ab (A-1-10) 
‘ Ka Ky a A. 


The solution of equation (A-1-10) for real values of h, and Hy implies 


L 
the restriction of the roots in a certain region. In general, if there 

are ' 4 '' independent feedback loops in a system, the coefficients of 

the characteristic equation consist of the terms hy» (hyhy),...¢h)..-h_). 
This leads to the following statement. 

If all the coefficients of the characteristic equation of a derivative 
controlled system are a linear function of the control parameter, then all 
feedback paths must be interconnected, (or touching). The above state- 
ment may be stated in another way: If all the coefficients of the character- 
istic equation of a derivative controlled system be a linear function of 
the control parameters, then those points from which the signals are 
taken as control variables cannot be used as feed-in nodes. 

Example A-1-4. In Fig. A-1-9, assume the signal flowing at the point 
marked "X'' can be measured, then it can be picked up as a control variable. 
Assume the feedback paths are as shown, they are interconnected. The 
characteristic equations are: 


Required characteristic equation = 5? + AS + ALS + A, = Q 


Characteristic equation A+b ab K KR 


a 
1 = h,S aks Kak, h, 
Ly = h, Ky (s°+ bs) 
The coefficient relations are: 
ponte 2 = elt by 
Kak, hy + eres — A,-ab 
Kak, = Ao 


h, and h, can be solved for arbitrary value of A's. 


1 2 
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(2) If zeros are involved in the forward path, the feedback 
from the intermediate nodes may introduce the problem of consistency. 
Consider Fig. A-1-10, the signal at "X" can be measured and then can be 


taken as control variable. The characteristic equations are: 


Required Chara, Equation oe + AS” + A,S + A, a0) 
Original loop 1 a+b abtk k, KK 
Lo ° hi S Kak, (S$ +o) kK, Ky h; Ka KX h, 
2 

L, 3 Kg (Sta) (S + bS)h, Km KaGstb)h, kab h, 
The coefficient relations are: 

dane" 

Ka K,h, + K, (4th, = A,- (atb 

akyh, + (4th, = A, (a+b) aoa, 


K,xb h, t oe K, xh, — A, —(abt Ka Ky) 
Ka Koc Tao 


Equation (A-1-12) are linear, but are not consistent for arbitrary 
A's, since the number of equations is one more than the number of unknowns 
(h, and hy). Substitute h, from the first equation into the second and 


third equation, obtain: 
KKh, t+ (CAtb) = A,-Catb) 
K,K,Xh, +X&b = A,— (ab+ kK Kp) 


Manipulate, get 


XA, -A, = AXA(KA+aA-3b) -Ab-k ky (A-1-13) 


Equation (A-1-13) is the restriction of the roots for this systm, i.e., if 
equation (A-1-12) has a solution, then equation (A-1-13) must be satisfied. 


Therefore, when a system has the required number of control parameters and 











Ka (Sta) 
GSD) 


Fig. A-1-11. Feedback of the output to more than one point. 


“ 
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the linear relationship of the coefficients, the choice of roots may still 
be subject to some restriction if the feedback is not proper. In general it 
can be shown that if the loop transmission has the same order(s) for both 
numerator and denominator, then it introduces the restriction of the roots. 
This leads to the following statement. 

| If the roots of derivative controlled systems are to be chosen arbi- 
trarily, then for each loop transmission, the order of the numerator must 

be at least one less than the denominator. 

(3) There is another obvious inconsistency which arises from 
feedback of the output to more than one point. Therefore, if an inter- 
mediate signal which is not followed by pure integrator, be fed to an- 
other point, the equations of the coefficients is inconsistent for arbi- 
trary choice of roots. Fig. A-1-11 is an example. 


The characteristic equation of Fig. A-1-1l is 


original loop 5° 57 5) ci 
, mo.) Qthtc ab+bC+ba = KK, R + abC 
ns Sa) CSc Ga 
i oo Kh, (ate)khy, bh, 
as 
ee ei NS K Khe 


There are three equations but two unknowns. This leads to a con- 
straint of the roots. 

(4) Since the feedback from an intermediate node is equivalent 
to derivative feedback from the output, then it is possible that the para- 
meters introduced by both control variables may have the same effect on 
the system; i.e., the parameters may not be independent. This can be 
easily detected by just looking at the loops. In-Fig. A-1-12, hy 
has the same effect as hy and hy» since hy is a combination of the first 


and second derivative of the output. That is, hy enters the coefficient 
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of the first and second derivative terms as h, and ho. Therefore in 
Fig. A-1-12, the system actually has only two control parameters. 
However, if the derivative of the signal at 'X'"' can be obtained 


and the scheme changes to Fig. A-1-13, then each parameter plays its own 


role, 
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Fig, A-1-12, Two parameters have the same effect, 





Fig. A-1-13. Each parameter has it's own effect 
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Die Summary. 

Derivative control is an effective method to control the coefficients 
of the characteristic equation of a system, since it uses derivatives of 
the system variables to control the coefficients of the characteristic 
equation. Since an intermediate signal feedback has the same effect as 
derivatives from the output, they are considered as in the category of 
derivatives. The criteria in Section I also applies to the intermediate 
signal feedback because of the condition (1) in Section 4. The inter- 
mediate signals do not introduce more control parameters but increase the 
possible ways of obtaining derivatives and the feedback path. 

In general, one constraint of the roots is almost always required in 
a system,that is, the product of the roots equal to the system forward 
gain, because the gain is usually specified from the stiffness of the 
system. This is only one degree of constraint, it does not introduce 
much difficulty. But if additional constraints are involved such as 
quadratic functions and inconsistent equations, then the solution is not 
so obvious. Therefore, it is logical to look for a set of linear and 
consistent control parameters. From the discussion in the previous sec- 
tions, it may be summarized as follows: 

(1) The number.of h parameters must be equal to (n-1l) 

(2) h - parameters must appear only in the coefficients from the 
first derivatives to the (n-1) derivative, 

(3) h-parameters must be linear combinations of the given co- 
efficients. 

Example A-1-5, Fig. 2-1-14 is a 6th order system. Assume two nodes 
between the energy-storage element are allowed to feed in or pickup signals, 


and the derivative of the signal at UX can be obtained. From the 
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' ; 6-2 7 
criteria of intermediate nodes, ( “A Oe 2) the system coefficients can 


be adjusted to arbitrary value by first and second derivative feedback, 
Fig. A-1-14 (a) (b) and (c) are the possible proper feedback paths. Fig. 


A-1-14(c) again has 5 possible combinations of H's. They are 


Hy ee h,S tales. Hy = eet Or hy) Hy = h,> aaa 
H, = h3S +h H, = h,S ths ane 
Hepaerie S14 nhs Ss” ) H, = hs +h.s H, = hs + hs 


Therefore, there are 8 possible ways to control the coefficients of the 
characteristic equation of this system, The choice of the scheme depends 


upon the specifications and perhaps the engineering judgments. 
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G, 









Fig. A-1-14. Example of A-1-5. 
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APPENDIX IT 
PARTITION OF A nTH ORDER SYSTEM WITH ONE ZERO FOR SINGLE SECTION CASCADE 
COMPENSATOR 


Take a 4th order system as shown in Fig. A-2-l. 






<(ew 


S(S*+Bs74 Bs +B.) 





Fig. A-2-l1. Cascade Compensator 


The controlled characteristic equation is 


S° + (B,+P)S* 4 (B+ Bp) S>+(@+PB)S° + BPs 
+ KAS + KX (2405 +Kpc =O 


S°+ (Bt p)St+ (Bt BPS + (B+ PB tka) Ss ate 
t(BP+KP +Kac)s +Kpc =O 


The characteristic equation of the uncompensated system is 
at > 2 ve 
Eres 4 (5,5 +(B+kK)JsS+ KC =O jhe Deny 


In order to express the coefficients of (A-2-1) in terms of the coefficients 


of the uncompensated system in a simpler form, define: 
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a 2 B, 


p, = 8. 
6 + K = B, (he= 3h) 
KC =. Bs 
Substitute (A-2-3) into (A-2-1) and manipulate, one obtains 
S° + (B,+p)S” + (B+ pB,) 5 + ce + Sou +t PB.) Ss” 
(A-2-4) 


i @eepmeoe,) Ser Bp == O 


Set the coefficients of (A-2-4) equal to the corresponding root-coeffici- 


ents and denote the arbitrary roots by < and Wr s one obtains 





© O O - | p 2(wW, - B; 
2 
B, O O | ae auor of LW, - Bs 
: ars aes iba 
B, Bey -2F Wn - Cy -B, + 22 | (A-2-5) 
2 
ee Bo -2fw, -W,, O e O 
2 
Ds O =O;. O O Sy Oo 
Where Co C, and C, are the coefficients of the reduced characteristic equa- 
tion. 


By Cramer's Rule. 


sap. Bo z z} 4 
F ? 3 2 Oy, 3 
[B- Bw, +2P Bw, + (4p) Jud” + 2s F296, + B,u,- Bw +2pu)” | 


Define: 


(A-2-6) 


Sr 

Y, se 
cs 
SA 
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a 


Te= L395 P| Ft ke 


For an uncompensated nth order system which has one open loop zero at 
(-C) and has the characteristic equation with coefficients B 1? B 9 

n- n- 
coe eB then 


= er 2 3 
BLE “9 Bln B Baty = Gy Bull,= ~~] ~ cL B4G% 8, tn * HPstdh ed 
2 5 2 
B, + a BW, y i 5, Ww, + sf Bn 4] + = (-$B,- 4.8, Un )EPur ty Bw | 


— 
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